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ABSTRACT 

Context. It is well known that the kinetic temperatures, T^, of the molecular clouds in the Galactic center region are higher than in 

typical disk clouds. However, the 7\i„ of the molecular complexes found at higher latitudes towards the giant molecular loops in the 

central region of the Galaxy is so far unknown. The gas of these high latitude molecular clouds (hereafter referred to as "halo clouds") 

is located in a region where the gas in the disk may interact with the gas in the halo in the Galactic center region. 

Aims. To derive in the molecular clouds at high latitude and understand the physical process responsible for the heating of the 

molecular gas both in the Central Molecular Zone (the concentration of molecular gas in the inner ~ 500 pc) and in the giant molecular 

loops. 

Methods. We measured the metastable inversion transitions of NH3 from (/, K) = (1, 1) to (6, 6) toward six positions selected through- 
out the Galactic central disk and halo. We used rotational diagrams and large velocity gradient (LVG) modeling to estimate the kinetic 
temperatures toward all the sources. We also observed other molecules like SiO, HNCO, CS, C 34 S, C is O, and 13 CO, to derive the 
densities and to trace different physical processes (shocks, photodissociation, dense gas) expected to dominate the heating of the 
molecular gas 

Results. We derive for the first time r kin of the high latitude clouds interacting with the disk in the Galactic center region. We find 
high rotational temperatures in all the observed positions. We derive two kinetic temperature components (~ 150 K and ~ 40 K) for 
the positions in the Central Molecular Zone, and only the warm kinetic temperature component for the clouds toward the giant molec- 
ular loops. The fractional abundances derived from the different molecules suggest that shocks provide the main heating mechanism 
throughout the Galactic center, also at high latitudes. 
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1. Introduction 

The interstellar molecular gas in the Galactic center (GC) region 
(i.e., in the inner ~ 1 kpc of the Galaxy) shows higher kinetic 
temperatures, 7\i n , than typical di sk clouds. Using m etastable 
inversion transitions of para-NH3, iGtisten et al.l (119811) derived 
kinetic temperatures in the range of 50-120 K towards Sgr A. 
I Mapping the (1,1), (2,2), and (3,3) inversion transitions of NH3, 
Mor ris et al.l 0983) found high kinetic temperatures (30-60 K) 
towards the denser portions of the GC region. Observing more 
high ly excited NH3 inv e rsion lin es, Mauersberg er et all ([l986a) 
and iHuttemei ster et al.l d!993bl) obtained kinetic temperatures 
7\j n > 100 K in all clouds in the GC including SgrB2 re- 
gion. Similarly high temperatures w ere also found in the cen - 
tral regions of nearby galaxies, (e.g., [Mauersber ger et al.l l2003). 
From metastable, i.e. J = K, inversion tran sitions of NH3 to- 
ward 3 6 clouds throughout the GC region, Hiittem eister et al.l 
(fl993ah suggested that in addition to a warm component there 
is also a "cool gas component" with 7\i n ~ 20 - 30 K. The ex- 
tended warm component in t he GC of ~ 200 K is not coupled 
with the dust (7d us t < 40 K. iRodriguez-Fernandez et aTl l2002; 
Odenwald & Fazio ll984tlCox & LaureiisffT9 89). High dust tem- 
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peratures (7d ust ~ 80 K) are only seen toward the SgrB2 molec- 
ular cloud, which is claimed t o be an anomalo us region, with 
recent massive star formation (Ball y et al.ll2010h . So far, to our 
knowledge, the T^ D of molecular clouds has never been deter- 
mined either at higher latitudes t owards the Giant Molecular 
Loops (GMLs, iFukui et aTl l2006h . or in the forbidden and/or 
high velocity components, explained by the barred potential 
model as XI orbits. 

The kinetic temperature of the molecular gas results from the 
balance of heating and cooling. Molecular clouds cool down by 
the collisional excitation of molecules and atoms followed by 
the ra diative emission of this energy from the cloud (Hollenbach 
1988). For the physical conditions pr esent in the GC, the cooling 
is dominated by H2 and CO, while Hiittemeist er et alJ (Q993a) 
propose that the dust in the GC region is also an important cool- 
ing agent. 

Dust heated via stars cannot heat gas sufficiently, just be- 
cause the gas is warmer than the dust. Several heating mecha- 
nisms for the G C region have been proposed, as e.g., heating 
by cosmic ra vs (IGtisten et alJ 1 198 lb iMorris et aTl [19 83): heat- 
ing by X-rays (Wat son et alFl 981; Nag avama et alj|2007l) : mag- 
netic ion-slipping ( Scald 1977 ). The dissipation of mechanical 
super sonic turbulence through shocks has been proposed for th e 
GC jFleckll981l;IWilson et alJ1982t lM"auersberg eret aljl9 86a). 
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The shocks can be produc ed by several phe nomena: supernova 
or hypernova explosions (Tanak aet al.ll2007|): resp onse of the 
gas in a barred potential model dBinnev et al.ll 1 99 lh ; and when 
the gas in the GMLs flows down their sides along the magnetic 
field lines, and joins with the gas layer of the Galactic plane gen- 
erating shock front at the "foot points" of the loops dFukui et al.l 
120061) . 

NH3 is one of the best thermometers for me asuring the gas 
kineti c temperatures in molecular clouds (see, iHo & Townesl 
119831) . Observing several metastable inversion transitions, one 
can determinate the kinetic temperature of the molecular clouds. 

In this paper, we derive for the first time the kinetic temper- 
atures of the molecular clouds in the disk-halo interaction re- 
gions (foot points of the GMLs and positions where the XI or- 
bits intercept X2 orbits in a barred potential). We use metastable 
inversion transitions of NH3 and other molecular tracers (SiO, 
HNCO, CS) to estimate the kinetic temperatures and densities, 
and discuss the heating mechanisms of the molecular gas in the 
GC. 



2. Observations 

2.1. Effelsberg observations 

We observed the metastable inversion transition of NH3 (J, K) — 
(1, 1), (2,2), (3,3), (4,4), (5,5), and (6,6) using the Effelsberg 
1 00m telescopeQ in April 20 1 0, and April 20 1 1 . We used the pri- 
mary focus A = 1.3 cm (18 - 26 GHz) receiver, which has 2 
linear polarizations, and a Fast Fourier Transform Spectrometer 
(FFTS) in the "broad IF band" mode with a bandwidth of 500 
MHz, providing an effective spectral resolution of 49.133 kHz 
or 0.386 km s -1 . We observed the (1, 1), (2,2) and (3,3) line 
simultaneously, with a band centered at 23.783 GHz, and the 
(4, 4) and (5, 5) line in a second setup (centered at 24.336 GHz). 
The (6,6) was observed in the third setup, centered at 25.056 
GHz, using the 100 MHz bandwidth FFTS spectrometer, which 
provides an effective spectral resolution of 9.827 kHz or 0.073 
km s _1 . The beam width of the telescope at 23.7 GHz is 42.2". 
The spectra were observed in position switching mode, with th e 
emission-free reference positions from Riquelme et al. (2010b), 
which were checked in the first setup, where the most intense 
lines are detected. Each position was observed for 12 min in the 
first setup, 24 min in the second setup, and 32 min in the third 
setup. The calibration in Effelsberg was done by the periodic in- 
jection of a constant signal (noise cal). To convert the data to 
T* K we corrected for the noise-cal (in K), opacity and elevation 
dependent antenna gairQ. The uncertainty in the calibration is of 
the order of 5 - 10%. The main beam temperatures, Tmb, were 
obtained by using Tmb = T* K ■ gij, where the beam efficiency, 
Beff, is 0.52 at 24 GHz. The pointing was checked every two 
hours against the source 1833-212, providing an accuracy of bet- 
ter than 10". 

In this work, we observe the positions selected in 
Riqu elme et al.l d2010al) . To avoid confusion, we use the no- 
tation of that wo r k. Th e "Central Molecular Zone" (CMZ, 
iMorris & SerabvrJ 1 19961) corresponds to the region about 
-0°.5 < / < 1°.5. Since the clouds of the CMZ are aligned along 
the Galactic plane within b~ 0°, this can be viewed as an exten- 
sion of the Galactic disk, towards galactocentric radii < 1 kpc 



1 Based on observations with the 100-m telescope of the Max- 
Planck-Institut fur Radioastronomie at Effelsberg 

2 http://www.mpifr-bonn.mpg.de/div/effelsberg/calibration/L3cmpf- 
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and will therefore be called "disk". When one observed position 
(from those called "disk") have kinematical components asso- 
ciated to both, the XI and the X2 orbits in the barred potential 
model, we called them explicitly as 'Disk XI" and "DiskX2". 
The source "Disk 2" which corresponds to SgrB2, is located to- 
ward the X2 orbits. Since this source does not have the velocity 
components associated to the XI orbits, we just call this source 
as "Disk". Gas toward the GMLs regions is labeled as "halo" 
in this paper, to differentiate them from the molecular clouds in 
the Galactic plane. This does not imply that the findings in this 
paper can be applied to the disk or the halo of the Galaxy as a 
whole, because all of the positions included in this work belong 
to the GC region. 

We observed six out of nine positions from Riqu elme et al.l 
d20T0ah visible from Effelsberg shown in Fig.Q] one in the foot- 
point of the GMLs (Halo 1), one in the top of the loop (Halo 4), 
two in the disk toward the location of the expected interactions 
between the XI and X2 orbits (Disk XI- 1, Disk XI -2, DiskX2- 
1, DiskX2-2,) in the barred potential model dBinnev et al.ll 19911) 
and a pair of positions toward the GC plane (Disk 1, Disk 2) as 
reference (Table [TJ. 

2.2. Observations with the IRAM 30m telescope 

In order to constrain the physical properties of the gas, we also 
observed the J - 2-1, y = rotational transitions of SiO, 
29 SiO, and 30 SiO, the J = 2- 1, 3 - 2 rotational transitions of CS 
and the J = 2 - 1 of C 34 S, the J = 10 - 9 transition of HNCO, 
and the J = 1-0 rotational transition of 13 CO and C 18 0. The 
observations were carried out with the IRAM-30m telescop^l 
at Pico Veleta (Spain) in several periods from June 2009 to 
October 2010. For the 3mm lines, we used the E090 band of 
the Eight Mixer Receiver (EMIRfl which provides a bandwidth 
of ~ 8 GHz simultaneously in both polarizations per sideband, 
and for CS (3 - 2) emission, we use the El 50 band of EMIR 
receiver, which provide a bandwidth of ~ 4 GHz simultaneously 
in both polarizations. As the backend, we used the Wideband 
Line Multiple Autocorrelator (WILMA), providing a resolution 
of 2 MHz or 6.6 km/s at 91 GHz and 4.1 km /s at 146 GHz 
We obs erved the nine selected positions from Riquelme et al.l 
(|20 1 Oal) that were all observable with the 30m telescope. In this 
work, we use the antenna temperature scale 7^, which can be 
converted to main beam temperature 7mb = T* k ■ ||^, where the 
forward efficiency Feff is 95% and the main-beam efficiency is 
Beff = 81% at 86 GHz, and Feff = 93% and Beff = 74% at 145 
GHz. The beam width of the telescope is 29" at 86 GHz, and 
16" at 145 GHz. 



3. Results 

Fig. [TJ shows the ammonia spectra taken toward each position 
in all the metastable inversion transitions observed in this work. 
Most of the metastable inversion transitions of NH3 were de- 
tected, except the (4,4), (5,5), and (6,6) of "Disk 1" and the (5,5) 
of "Halo 4". The criteria used to define if a emission line is de- 
tected or not, was to have a line peak temperature > 3cr rms , where 
is the root mean square per spectral channel. If the intensity 
of the line does not reach this value, we still assume that a line 

3 Based on observations carried out with the IRAM 30m telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany), and 
IGN (Spain). 

4 http://www.iram.es/IRAMES/mainWiki/EmirforAstronomers 
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Table 1. Observed positions in NH3 lines. 
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Fig. 1. Spectra toward selected positions in the GC in the metastable inversion tran sitions from (1,1) to (6,6 ) of NH3. The pos itions 
are indicated in the HCO + integrated intensity map from iRiquelme et aU d2010bl) . The GMLs found by iFukui et alJ d2006l) are 
indicated in blue. The positions which could not seen from Effelsberg are indicated with green crosses. As indicated in TableQ] our 
Disk 2 position corresponds to SgrB2. 



is actually detected if the line has an integrated intensity in the 
velocity width (as defined by the (3,3) line which presents the 
highest signal-to-noise ratio) > 3cr. 



3.1. Optical depth of NH 3 

Each NH3 inversion transition is split into five components: a 
"main component" and four symmetrically placed "satellites" 
(the quadrupole hyperfine (HF) structure). Due to the large 
linewidth of the molecular clouds in the GC, the magnetic split- 
ting (< 0.2 km s _I ) cannot be resolved. Under the assumption of 
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local thermodynamical equilibrium (LTE), the relative intensi- 
ties of the four satellite HF components can be used to estimate 
the optical depth t of the main component of the metastable in- 
version transitions. Knowing t, we can estimate the NH3 col- 
umn density, and the rotational temperature from the ratios of 
the peak or integrated intensities. 

We use the "NH 3 method" from CLASS0 to determine the 
optical depth for the (1, 1), (2, 2) and (3, 3) lines. To define the 
linewidth (which was used as a fixed parameter in the NH3 
method), we use the (3,3) transition, because these spectra 
have the best signal-to-noise ratio in our observations and the 
HF components are much weaker than those of the (1,1) and 
(2, 2) lines. As we can see in Table IA.4I all the NH3 lines ob- 
served in this work are optically thin toward all sources, ex- 
cept the (1,1) trans i tion tow ard SgrB2. Following the criteria of 
Hiittemeis ter et alJ dl993ah based on the lower peak intensities 
in these lines, with respect to the (1, 1), (2, 2) and (3, 3) ones, we 
assume that the (4, 4) and (5, 5) are also optically thin. Table lA.4l 
presents the results from simple Gaussian fits for all the observed 
positions, allowing all the parameters to be free. 

3.2. Physical conditions of the gas from CS and A/H 3 

To derive the nQiz) and 7\i n of the gas, we combine the 
CS and NH3 molecular emission, in an iterative way. First, 
we used MASS A softwar^l to derive the rotational temper- 
atures and column dens ities using Boltzmann diagrams (see, 
iGoldsmith & Langell999] for a detailed explanation and equa- 
tions of the method) (Tabled- The rotational temperature, which 
is a lower limit of the actual kinetic temperature, T^ B , was u sed 
as a fixed parameter in RADEX (see, Ivan der Tak et alj2007l for 
detailed explanation of the formalism adopted in this statistical 
equilibrium radiative transfer code) to derive the n(Hi) and CS 
column densities. Then, using the n(H2) obtained from CS, we 
used RADEX to derive the kinetic temperature from the para- 
NH3 transitions (see 33.2.21 ). With the kinetic temperature, we 
derived then the final n(Yii) and CS column densities (Table[3]l. 

3.2.1 . n(H 2 ) derived from the CS data 

We used the non-LTE excitation radiative transfer code RADEX 
to derive the «(H2) and CS column densities from line intensi- 
ties of the observed CS lines. The modeling suggests that the 
CS emission is optically thin with opacities ranging from 0.05 
to 0.96. The results are showed in Table [3] The error were esti- 
mated assuming a 10% calibration error as the typical flux cal- 
ibration uncertainty at the 30-m telescope, and we give an up- 
per and lower value based on the minimum and maximum value 
from the LVG diagrams (see from Fig. |B.9| to Fig. |B.19l ). It is im- 
portant to note that n(Hz) in some sources is poorly constrained, 
which translates in large errors or upper limits shown in the Table 
|U If we derive the n(H2) using the rotational temperature (which 
is a lower limit to the kinetic temperature), the n(H2) differ on 
average by ~ 27%. 



and the velocity widths (see Table |A~41 >. we can derive the 7\j n 
and /Vnh,- Fig. |2] shows an example of this procedure, and Table 
|4]shows the results. In Fig.|2]and from Fig. IB.ll to Fig. IB. 81 we 
show in blue the results corresponding to the metastable inver- 
sion transitions (1, l)-(2,2) (low rotational temperature), and in 
red, the results corresponding to the metastable inversion tran- 
sitions (2,2)-(4,4)-(5,5) (high rotational temperature). For the 
cases where only one temperature regimen was a possible so- 
lution, we plotted the result in red in the LVG plot. LVG mod- 
els indicate that the results from LTE are reliable. Additionally, 
for every observed position, we checked the two temperature 
component assumption by comparison to synthetic spectra un- 
der LTE approach using MASSA software. We found that for 
the positions where we derived two kinetic temperature compo- 
nents, the modeled line profile fit better the observed emission, 
while a single warm component was not enough to reproduce 
the observed profile. When the (4,4) or (5,5) inversion transi- 
tions were not detected, the upper limits to their emission were 
plotted in dashed lines. This upper limit to the emission was ob- 
tained as 3o" lms level. The individual fits to all sources are shown 
in the Online Appendix (Fig lB.ll to lB.81 l. 

As a result of our analysis, we derive two kinetic tempera- 
ture (one cool and one warm) components in the CMZ and only 
one warm component in the halo positions. In the CMZ the cool 
component range from 23 K to 68 K with an average value of 
40 K; and the warm component range from 80 K to > 300 K 
with an average value of 150 K. This reference values should 
be taken with caution due to the large uncertainty of the kinetic 
temperature derived from the LVG (see below). To estimate the 
uncertainty of the derived parameters, we computed the^ 2 of the 
line intensities over the grid used for the LVG model. We impose 
Ay 2 = y 2 - v 2 . = 1, which translates in the 68.3% confidence 
level projected for each parameter axis (see, e.g., iPress et al.l 
1992, Section §15.6). The black ellipses shows the error in the 
model (Fig. 



3.3. Column densities and relative abundances from other 
molecules 

To shed light on the physical processes that are heating the 
molecular gas, we derived also relative abundances of NH 3 
with respect to the f ollowing molecules: SiO, w hich is a well- 
known shock tracer ( Martfn-Pinta do et aD 19921). CS which is a 



high-density gas tracer (n > 1 4 cm 3 Mau ersberge r & Henkell 
1989; iMauersberger et all 1 9891) . and HNCO whic h is a tracer of 
shocks, and very high densities, (n^ > 10 6 cm~ 3 I Jackson et al.l 
Il984t iMartm et al.l 120081: Hinchenko etaD l2000h with a high 
photodissociation rate (Table lATTT ). We also derive the fractional 
abundances of these molecules with respect to H2 as traced by 
13 CO (Table|6](. We assumed that they arise from the same vol- 
ume. This assumption may not be fulfilled for all of our observed 
position due to the observed differences in the velocity center 
and linewidth (see Table lA.5l ). In all the calculations we assume 



that the GC sources are extended, therefore we take ~ Tmb ■ 



3.2.2. LVG analysis from NH 3 

To estimate the kinetic temperatures of the gas, we also used a 
non-LTE excitation and radiative transfer code RADEX. Using 
the value of n(H2) derived from the CS LVG analysis (Table @J, 

5 http://www.iram.fr/IRAMFR/GILDAS 

6 http://damir.iem.csic.es/mediawiki- 1 . 12.0/index.php/MASSA 
_Users_Manual 



3.3.1 . Column densities of SiO, HNCO, 13 CO, and C 18 

We used the non-LTE excitation radiative transfer code RADEX 
to derive the column densities. For the species with only one ob- 
served transition (SiO, HNCO, and CO isotopes), we are forced 
to make some assumptions about the physical properties of the 
gas (7\j n and n(H2)). We used the kinetic temperatures derived 
by NH 3 , and the n(H 2 ) from the CS data (Table 0). The re- 
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Table 2. NH3 physical parameters (rotational temperatures and column densities) derived for each source using MASSA software. 
Bold faced values indicated the most likely result consistent with the non-LTE analysis. 



Source n" T m T mt T ml T Iot /V(NH 3 ) 2V(NH 3 ) /V(NH 3 ) yV(o-NH 3 ) 



(11-22) (22-44-55) (11-22-44-55) (33-66) (11-22) (22-44-55) (11-22-44-55) (33-66) 







[K] 




[K] 




[K] 




[K] 


10 14 cm" 2 


10 14 cm" 2 


10 14 cm" 2 


10 14 cm" 2 


Halol 


1 


46.7 ± 


0.9 


117.5 i 


1.7 


92.7 ± 


0.8 


81.0 + 0.8 


2.14 ±0.07 


3.53 ±0.10 


3.72 ± 0.06 


5.37 ± 1.1 




2 


43.2 ± 


0.3 


131.8 ± 


0.9 


96.6 ± 


0.3 


156.4 ±9.6 


4.93 ± 0.06 


9.09 ±0.11 


9.49 ± 0.06 


11.1 ± 1.2 




3 










80.9 ± 


1.1 


169 ± 31 






2.06 ± 0.06 


3.3 ± 1.0 


Halo4 












56.5 ± 


2.3 


138 ± 10 






5.48 ±0.45 


5.39 ±0.76 


DiskXl-1 




27.2 ± 


1.5 


156 ± 


12 






164.1 ±6.6 


4.85 ±0.53 


7.7 ± 1.2 




12.07 ± 0.91 


DiskX2-l 




27.4 ± 


0.8 


91.7 ± 


3.6 






122.8 ± 3.8 


9.75 ± 0.52 


9.87 ± 0.78 




16.56 ± 0.98 


DiskXl-2 




33.1 ± 


2.2 


112.7 i 


6.5 






138.2 ±5.1 


3.99 ± 0.48 


5.75 ±0.71 




10.93 ± 0.80 


DiskX2-2 




22.0 ± 


1.4 


72.5 ± 


6.1 






106.9 ±9.7 


4.44 ± 0.55 


2.85 ±0.58 




5.09 ± 0.90 


Diskl 




19.1 ± 


2.9 


154 ± 


48 






139 ±26 


0.74 ± 0.20 


0.66 ±0.41 




0.93 ± 0.39 




2 


26.9 ± 


2.4 


82.9 ± 


9.8 






100 ± 12 


1.41 ±0.20 


1.29 ±0.34 




2.02 ± 0.46 


Disk2 


1 


38.1 ± 


0.6 


110.3 4 


1.0 


90.6 ± 


0.5 


95.7 ± 1.5 


31.6 ±1.0 


48.1 ± 1.1 


60.3 ± 1.0 


61.0 ±2.0 




2 


50.8 ± 


1.6 


145.6 4 


2.3 


113.8 4 


0.5 


112.9 ±0.6 


42.1 ±2.1 


84.9 ±3.1 


101.6 ± 1.0 


133.2 ± 1.4 




3 


33.0 ± 


0.9 


63.1 ± 


1.2 


51.7 ± 


0.5 


52.4 ± 3.2 


33.1 ± 1.3 


32.0 ± 1.7 


47.2 ± 1.0 


37.6 ± 6.4 



" Cloud number denned by the different velocity components (see Table lA.4l 

Nipara - NH 3 ) correspond to the sum of all observed para-NH 3 column densities. When the data is consistent with a two temperature model, 
N(para - NH 3 ) correspond to the sum of column 7 and 8. If only one temperature regime is present, the N(para - NH 3 ) corresponds to column 
9. 



Table 3. Physical parameters derived from CS using non-LTE (RADEX) model. The kinetic temperature used in RADEX is indi- 
cated. When more than one T^ m regime was present in one position, we derived the physical parameters from each T^\ n . We assumed 
a 10% calibration error. 



Source 



Cloud 
number 



Tkin 

[K] 



Tex CS(2-1) 
[K] 



Tex CS(3-2) 
[K] 



T CS (2-1) Tcs (3-2) 



Halol 



115 

90 
135 



< 3.8 
12 6 93 
4 6' 7 



3.8X1 



< 3.6 
8.2* 



>0.25 

25 35 

015 

31'* 



22 u:vi 
0.28'" 



>0.16 

54° 34 
35 077 



0.15=| 



Halo2 



113 
113 



4.0 
3.8° 



Halo3 
Halo4 



113 
113 

95 



4.7: 

6 l 3 ' 



4 4 1J 

5 2 L2 
4.7il 



0-61^' 
07°" 

• Q.Q4 

0.21° 28 



0.75 
0.05° 

0.17° 



Halo5 



95 



10.4? 



7. Of 



O.m 



0.21" 



Disk-Xl-1 



38 
300 



10.9- 
9J& 



"•'1.6 

5 9 2 5 
J - y \ 4 



7.1 
6.8 



1.4 



5 311 

5.o!- 3 



"■ i -'0.0') 



"037031- 

0.50' 23 



0.24 



25°'^ 



-OA^r 

0.55 ; 85 



Disk-X2-1 



38 
100 



13 W~ 

12.9! 5 , 5 



0.43" 23 



Disk-Xl-2 



52 
215 



8.4 



0-16; 

0.18°?! 



Disk-X2-2 



28 
95 



lO^T 

9 3" 



0.20 
0.22 



0.14 



0.30; 

0.31 



S9 
0.11 



Diskl 



23 
154 

38 
82 



< 3.4 

< 3.0 
5 5 2 ' 
5.0 : ? 



< 3.3 

< 3.0 
4 7'° 



4.5 



o.x 



5 3° 

5.7:1 

9.4] 

9 3 3:s 



>0.22 >0.14 



>0.57 
0.17" 20 
0.21° 



§7 

1 0.12 



16*™ 
20° 28 

70' 4 - 8 
'"0.41 

19° 22 
19° 27 



>0.30 
17°'° 

0.05 

19° 18 



21 UI2 
23 017 
96'" 
45° 23 
45° 26 



Disk2 



68 
200 
145 
50 
80 



u,1 1.6 

6.7- 
16.1| 
16.11" 



suits are shown in Table [5] The column densities agree within 
a factor of 3-4 if we use the LTE approach (Table IA.5b for a 
T ex = 10 K (which i s con sistent with the value derived for SiO by 
Hiittemeister et al. ( 1998), and with our T ex value derived from 
CS shown in Table[3]l. 



Using n(H2) ~ 10 3 cm" 3 (because the critical density of CO 
is lower than for CS) the column densities for the CO isotopes 
are a factor 2-4 lower than using the n(H2) ~ 10 4 cm -3 . 

The total column density of H2 can be estimated from the 
CO isotopologues, that have lower optical depths than the main 
isotope, and therefore a more reliable estimation of the col- 
umn density (iV(H 2 ) = N( x C y O) x [ x C y O/H 2 ], where x and 
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Table 4. Kinetic temperatures and densities derived from LVG calculations from NH3 and CS data 



Source 


Cloud 
number 


low temperature 

Zkin «(H 2 ) 


high temperature 
Tkin «(H 2 ) 


single temperature from p-NH? 

Tkin «( H 2) 






r \r 1 


1 n4 „ m — 3 

10 cm 




1U cm 




1U cm 


naio 1 


1 
1 

2 

3 










> Li J 

v. on 
1 

> ijj 


< l.UU 

/iq5.63 

i.DS 123 


Halo2 


1 

2 










1 1 ic 


1 26 1 - 56 
i./b 098 

^ o ^ 1 


Halo3 


1 

2 










1 1 J 

113'' 


1 50 1U5 
3 55 2 - 76 




Halo4 


1 










95 


^•' J 1.49 


Halo5 


1 










9S r 


7 47 5 -" 


DiskXl-1 


1 


38 


13 s " 1 


> 300" 








DiskX2-l 


1 


38 




100 








DiskXl-2 


1 


52 


13.3^(4 


215 


|g| 






DiskX2-2 


1 


28 


14 |H 


95 








Diskl 


1 


23 


< 1.12 


> 154" 


< 0.16 








2 


38 


4.57 3 - 84 


> 82 6 


251 2.40 






Disk2 


1 

2 
3 


68 
50 


4. O^-iS 
^■ ?/ 2.40 

15 ' 8 M4 


200 
80 


7 472.WI 

11 6 8 - 32 


> 145 


2 5 J2.50 



" LVG gives a 7\in greater than 300 K which is the value allowed by the collisional rates given bv lDanbvetaIl in"988). 

* Due that the values for the (4,4)~(5,5) metastable inversion transitions are upper limits to the actual T MB , the modeled curve in the LVG plot 

is outside the allowed range, and we give a lower limit to the kinetic temperature using the rotational temperature form the LTE plot. 

c The assumed value for the kinetic temperature for positions were there are no NH 3 observations is the average of the T^ n from similar 

positions. 



Table 5. Column densities from different molecules derived from RADEX using as a fixed parameter the kinetic temperatures and 
the n(H 2 ) (Table H 



Source 



Cloud 
number 



Tldn 

[KJ 



N(p-NH 3 ) 
[10 I4 cm" 2 ] 
— 1 07 im 

A - v ' n 27 
7 740:11 
^'^026 

70° 15 



N(CS) 
[10 13 cnT 2 ] 



> 6.73 
15.99|» 
6 53 17 - 17 

c 7s: 16.65 

2.73 

> 2.91 



N(SiO) 
[10 13 cm- 2 ] 



N(HNCO) 
[10 13 cm" 2 ] 



N(C 34 S) 
[10 I3 cm" 2 ] 



N( 1J CO) 
[10 16 cm- 2 ] 
4.47^ 

4 03^ 
2 00 - 75 



N(C ls O) 
[10 16 cnr 2 ] 
— 012™ — 
97^ : " 
16 008 



N(H 2 ) fl 
[10 21 cnr 2 ] 
23 53 r ' 4 
21 21 353 

3 95 



Halol 



115 

90 
135 



> 2.24 
4 35 1% 
1.88*2 



>0.50 
2 89 - 79 



1 20 044 



1 64 4li3 



33 l b/ 

"■ JJ 0.21 

> 0.13 



Halo2 



113 
113 



"^2 _ 52° 5_ 

J - 71 1.07 

12 64 3 ;i3 
Ml 

13.78 ?s? 



76 2 35 

/u 0.38 

>0.34 



"198^ 



H - H/ 1.74 

1 27 - 72 



1 43<n5 

29 034 
——242- 
40" 42 

0.47"-' 4 



Halo3 



113 
113 



Too™" 



Tog 35 " 



735™ 

2.01 



25 U - Uii 



38.68^ 



Halo4 



95 



1 99° 



U 0.38 

0.99^ 42 



7 coo:v6 

4o%r 



Halo5 



95 



L99JI 

^'^'0.69 

3 30 043 
1 90" 51 

1 00° 60 

m 

40 1 84 

"' w 026 

>0.16 

63° 20 
>0.45 



"■"in 
3 91™ 

9- 



64" ' 

fH? 

°- 71 18 



8.84|;| 2 

32.8o 75 

16 9^ 

29 2 7 - 55 

"IToTfT 

T2-sy^ 

28.7"* 



51 u - lg 

Q-18 

Lij7 0S7 

1 37033 

J 0.27 

2 51°- 77 
-032^- 

1 14W3 

2 62"- 52 

1 

19™ 



46 53 5 ^ 
172 63 4 * 42 
33.^ 

42.58^ 3 

57.40!«i 2 
25 2i ' i)f 

Z64*! 
4 52* ):7 ° 



Disk-Xl-1 



38 
300 



3.98, 



19.07'^ 336 



Q ^qioo 
7 ""l.02 

10 6 130 

2.01 



74 22 
"■'^0.11 

93" * 
"■"0.15 

96" 48 
>0.03 
>0.09 
16 013 
\9 m 

"ToT^ 

J -^ J 1.04 

> 5.67 



Disk-X2-1 



38 
100 



"• J/ 5.44 

22.56 37 f 
10.60-- 

10.67|* 
1 1 07 s:2S 

> 2.90 

> 7.20 
3.89 280 



'1.12 

3.60 676 

2 10°- 78 

"122™" 

27™ 
"■ z/ oo8 

>0.25 

> 1.43 
0.45 037 
50^ : ^ 

316^ 

j. iu 106 

> 5.63 



o iq0.8i) 
2 96 % 

_P39 



Disk-Xl-2 



52 
215 



2.00^_ 



Disk-X2-2 



28 
95 



> 1.12 

> 2.82 
2 72°- 79 

2!oo! : i 



1 26 i):§2 
Q-55 

1 32°-" 

2 26 8: ^ 
3-46?l 



Diskl 



23 
154 

38 
82 



4.54 



TO 7 ™" 
14.1™ 



> 16.51 
16.89t; 3 ^ 



Disk2 



68 
200 
145 
50 
80 



25.1 
14.1 



> 264.4 



101.28 2847 



19 63 648 



1 0.05 

12 l 07 
iZ " 1 007 



17.05S; 



20 2 56 56 94 

156.60 



" Column densities of H 2 derived from 13 CO using a conversion factor of 5.0 x 10 6 {Rodriguez-Fernandez e t al.l200lh for the normal GC gas, 
and a factor of 1.9 x 10~ 6 for the gas in the disk-halo interaction regions (see text for details). 
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Radex, para-NH 3 Halo 4 Radex, para-NH 3 Disk XI -2 



13 



12 




NH,(5,5) 
NH,(4,4) - 



n(H 3 ) = 2.75xl0 4 cm" 
AX 2 = 1 



16 r 



f)0 



100 



150 
T»«. [K] 



13.0 



11.0 









Halo 


4 


































N 







200 

Eu/k (K) 




13 - 



12 t 



250 300 



5!) 



n(H 2 ) = 12.59x10' cm" 
n(H 2 )=5. 62x10" cm" 3 
AX 2 = 1 



150 
T kl „ [K] 



13.0 



' 12.0 









-2 

































200 

Eu/k (K) 



Fig. 2. Example of rotational (bottom) and LVG (top) diagrams of NH3. Left: Halo4. Right: DiskXl-2. In the Boltzmann diagram 
(bottom), para-NH^ species are fitted with the red line, and the ortho-NH3 species are fitted with the blue line. In the LVG diagrams 
(top), we plot the peak intensity of the metastable inversion transitions of para-NH^. For the source Halo 4, we derive only one 
kinetic temperature component (warm), which is plotted in red lines, and the error in red shaded region. The source Disk XI -2 
was modeled using two kinetic temperature component model, one warm and one cold, that were plotted in red and blue lines 
respectively with the correspondly associated error shown as shaded regions. We show the n(H2) used in the LVG model for the 
warm and cold components. To estimate the uncertainty of the derived parameters, we computed the^ 2 of the line intensities over 
the grid used for the LVG model. We impose Ay 2 = x 2 _ ^min = ^> w hich translates in the 68.3% (1 <x) confidence level projected 
for each parameter axis 
(see text for details). 



y correspond to the isotopic substitution used for the carbon 
and oxygen atoms). For our calcul ations, we assume an abun- 
dance ratio CO/H2 of 10 (see, e.g., iRodrfguez-Fernandez et al. 
I2001L and references therein), and we use the 13 CO emission. 
Therefore, we a lso nee d the 12 C/ 13 C isotopic ratio. Recently, 
Riqu elme et alJ (1201 Oah derived a high 12 C/ 13 C isotopic value 
(> 40) in some of the sources studied in this work. Such a high 
isotopic value was found mainly toward the disk-halo interac- 
tion r egions, therefore, we still u se the standard value of 20 (see, 
e.g., I Wilson & Matteucci|[l992l) in the typical GC gas for the 
sources "Diskl", "Disk2" and for the sources with kinematic of 
X2 orbits (DiskX2-l, DiskX2-2). For the sources which are in 
the disk-halo interaction regions, we used the value of 53 cor- 
respo nding to the typical v alue found in the 4 kpc molecular 
ring (Wils on & Rood! 19941). w hich also was used by Tor ii et al.l 
d2010al) and lKudo et all d201 ll) in the GMLs regions. This trans- 
lates in using a [ 13 CO/H2] conversion factor of 5.0 x 10" 6 for 
the normal GC gas, and 1.9 x 10~ 6 for the gas in the disk-halo 
interaction regions (Table |5). We decided not to use the C 18 
emission to trace the total column density of H2, because of the 
uncertainties in the 16 0/ 18 isotopic ratio in the disk-halo in- 
teraction regions, which could be affected by the unprocessed 



gas th at is being accreted toward the GC region ( Riquelm e et al.l 
I2010al) . 



3.3.2. Fractional abundances 

We derived beam averaged fractional abundances with respect 
to H2 for all the observed molecules (Table and in Table lA.il 
we show the results of the fractional abundances of NH3 with 
respect to the other molecules, and the fractional abundances of 
SiO and HNCO with respect to CS and C 34 S. 

The results for Disk 1, will be disscused below. X(SiO) varies 
from 0.07 - 4.17 x 10~ 9 with the largest value toward the foot 
points of the GMLs. X(HNCO) shows less variation in all the 
observed sources ranging from 0.28 to 3.08 x 10" 9 , except for 
the Disk 2, which has a large abundance of 17.31 x 10~ 9 . The 
fractional abundances with respect to H2 depend on a reliable 
estimation of the H2 column density, which depends in a num- 
ber of assumption (H2 to CO conversion factor, physical param- 
eters used to derive the column densities of 13 CO). Therefore, 
we also compare the column density of the different molecu- 
lar tracer (shock, photodissociation) with respect to CS (dense 
tracer) and its C 34 S isotope (Table I A. It . which is almost cer- 
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Table 6. Fractional abundances of SiO, HNCO, NH 3 , CS, and C 34 S with respect to H 2 



Source 


n" 


X(SiO) 


X(HNCO) 


X(NH 3 ) 6 


Z(CS) 


X(C j4 S) 






xl0~ 9 


xlO -9 


xlO s 


xlO" 9 


xlO 10 


Halol 


1 


1.42 ± 0.95 


0.32 ± 0.21 


1.70 ± 1.25 


4.27 ± 2.84 






2 


2.33 ± 0.79 


1.45 ± 0.39 


2.12 ± 0.40 


9.16 ± 2.25 


6.36 ± 1.84 




3 


4.17 ± 3.73 




1.42 ± 0.77 


1 A t C , 1 O O O 

14.15 ± 12.38 




Halo2 


1 

ri 

I 












Haloj 


1 

L 












xiaio4 


i 
t 


U.j4 ± U.ZU 


Y.Ly + U.jo 


i i i i n 
1.11+ U.J / 


i n/i _i_ 1 1 3 

Z.U4 + 1. 1 J 


L.J 1 ± U.o / 


Halo5 


1 












DiskXl-1 


1 


0.39 ±0.12 


0.28 ± 0.07 


0.35 ± 0.09 


1.31 ±0.37 


0.67 ± 0.20 


UlSKAz-i 


t 
i 


1.2U ± U.JB 


i no i 1 1 o 
3. Do ± 1.12 


1. /Z ± 0.41 


6.94 ± J. 24 


2.^4 ± l.Jy 


DiskXl-2 


1 


0.77 ±0.18 


1.12 ±0.25 


1.10 ± 0.21 


3.95 ±0.86 


2.60 ± 0.54 


DiskX2-2 


1 


0.07 ± 0.02 


0.52 ±0.10 


0.76 ±0.18 


2.91 ±0.70 


2.55 ± 0.60 


Diskl 


1 


3.60 ± 1.01 


8.44 ± 2.37 


2.40 ± 0.67 


21.63 ±6.08 


2.57 ± 0.72 




2 


1.12 ±0.49 


4.22 ± 1.45 


1.94 ±0.36 


9.94 ±4.13 


4.21 ± 1.87 


Disk2 


1 














2 
3 


0.37 ± 0.26 


17.31 ± 12.10 


3.29 ±2.30 


1.08 ±0.76 


3.71 ±2.60 



" Cloud number denned by the different velocity components. 

b Column density of NH 3 correspond to the N(p - NH 3 ) + N(o - NH 3 ) = 2 x N(p - NH 3 ) assuming a ortho-to-para ratio of 1. 



tainly optically thin becaus e in the GC regio n 3 S is ~ 10 times 
rarer than the main isotope (IChin et alj|1996l) . 

Although there are no big differences from source to source 
in the jV(SiO)/W(CS) and N(SiO)/7V(C 34 S) ratios, we can see 
that the highest values are found toward Halol and Disk XI 
sources, with a difference up to 1 order of magnitude if we 
compare the DiskXl-2 with the DiskX2-2 sources. The rela- 
tive abundance jV(HNCO)/jV(CS) and jV(HNCO)/7V(C 34 S) to- 
ward the Disk 2 source is by far high est. The source Disk 2 corre- 
sponds to SgrB2M (20",100") from lMartm et alJ (120081) . which 
is classified as a "typical Galactic center cloud" in their work. 
They found a large HNCO/ 13 CS abundance ratio of 68 ± 13 in 
that source. 

We exclude the source Disk 1 from the previous analysis, 
because the determination of their physical parameters (kinetic 
temperature) would be overestimated. The metastable inversion 
transitions (4,4) and (5,5) of NH3 were not detected, therefore 
the kinetic temperature should not be high. In our radiative trans- 
fer calculations we use an upper limit to the rotational tempera- 
ture (154 K), which was taken as the kinetic temperature of the 
gas. It is probable that the actual kinetic temperature (if there is 
a high temperature regimen in this source) could be much lower 
(similar to the value for the disk X2 positions, which correspond 
to typical gas in the CMZ). The LVG column density of SiO on 
this source (Table|5]) is a factor of ~ 30 larger than the value from 
the LTE (Table lA~5l l. while the differences for the other positions 
are only a factor 2-4. 



4. Discussion 

4.1. Kinetic temperatures toward the Galactic center region 

The derived kinetic temperatures for the halo positions are con- 
sistent with one high temperature regime (> 90 K), while the 
clouds in the CMZ are consistent with two temperature regimes 
(~ 40 and ~ 150 K). 



4.1 .1 . Single temperature regime in the loop region 

We derived a single high kinetic temperature regime (> 95 K) 
for the hal o sources (Halo 1 a nd Halo 4), towards the GML dis- 
covered bv lFukui et alJ ([2006). Surprisingly, both the Halo 1 po- 
sition in the footpoint of the GMLs and Halo 4 in the top of the 
loop do not show any trace of low kinetic temperatures, which 
otherwise a re present throughout the CMZ a s discussed in previ- 
ous works (iHiittemeister et al.ll 1 993aL 1 1 998h . iTorii et all (l2010ah 
derived kinetic temperature of 30-100 K or higher, and den- 
sities of 10 3 ctrr 3 using multitransitional CO observations to- 
ward the foot point of the GML (loop 1 and 2). This foot point 
corresponds to our Halo 2, and Halo 3, which could not be ob- 
served with t he Effelsber g telesco pe due to their low declination. 
Furthermore, iTorii et all (1201 Obi) made a comparison between 
the foot points 1 and 2 with the two broad velocity features, the 
Clump 2 and / = 5°. 5 complex, finding that they share common 
properties such as the vertical elongation to the plane and large 
velocity spans of 50 - 150 kms -1 suggesting that they have a 
similar origin. Therefore, the physical processes that are occur- 
ring in the Halo 1 position should be similar to those in the well 
studied foot point of the loopl and 2. 

Is our Halo 1 position really located at the foot point of the 
GML?, or it is just along the line of sight given that we see the 
GC region edge-on? This source is at I = 5°. 5 and at ~ 790 pc in 
projection of the GC ( /, b) = (0, 0) positio n assuming a GC dis- 
tance of 8.23 +0.2 kpc (Genz elet alJ2010[). It was previously ob- 
served by several author s (Bitran et al. 1997; Fukui et al. 2006; 
ISawada et alJ200lHLee & Leel2003HRiquelme et al.ll2010bl) : its 
main velocity component is at 98 km s _1 (from the HCO + data 
from Riquelm e et alj2010bl) . and a large velocity width in all the 
species is observed in this and previous works, which indicates 
that this source indeed belongs to the GC region. Furthermore, 
we cannot rule out that the gas seen in the Halo 1 position has 
some velocity component belonging to the GC region, but at 
smaller or larger distances. High resolution maps of the foot 
point region as well as maps of the complete loop are needed 
to reveal the morphology and kinematic of the complete loop to 
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confirm the association of this position to t he GMLs s c enario , 
because this loop is tentatively detected by iFukui et all d2 006). 
Therefo re, the high 12 / 13 C isotopic ratio found by Riquelme et al. 
(2010a) toward this position, and towards the well studied foot 
point of the loop 1 and 2 would provide evidence for the GMLs 
scenario. 

Additional support for the GMLs scenario comes from 
the kinetic temperature gradient and large NH3 abundances 
of the high metastable inversion transitions. The sense of the 
temperature gradient can help to establish whether the shocks 
are due to the GMLs scenario or to the ejection of gas from 
the disk due to star formation. Temperature that decreases 
from the disk (low latitudes) to the halo in the GC region 
would indicate that the material is falling from the halo to 
the galactic disk supporting the GML scenario, because the 
post shocked gas which has cooled down is at larger latitude 
than the re cently heated material at the shock front (see e.g., 
Fig. 16 of lGenzell[l992l) . The gradient will be in the opposite 
way if the material is being ejected. We observed that the 
kinetic temperature is slightly larger in the foot point than in 
the top of the loop, which tentatively supports the loops scenario. 



4.1 .2. The two temperature components model in the CMZ 
clouds 

In the CMZ (Disk XI- 1, Disk XI -2, DiskX2-l, DiskX2-2, and 
Disk 2), our results are consistent with a two components model, 
with at least at two different temperatures, one cool and one 
warm (see Table 2] and ffl- Thi s result is in agreement with 
IHiittemeister et al.l7l993aTTl998l) . who found that the data were 
consistent with two rotational temperature components: one cool 
(~ 25 K) and other warm (> 100 K). Furthermore, they found 
that for a typical GC molecular cloud, 25% of the gas has high 
temperatures, and this gas has low H2 density; while the remain- 
ing 75% of the total gas mass is cooler at densities of > 10 4 



cm . Both gas components are in pressure equilibrium. Our re- 
sult, on the other hand, indicates that there is as much gas in 
the low temperature component (~ 50%) as in the high tem- 
perature component (~ 50%). A possible explanation for such 
a discrepancy may be a selection effect: The positions from 
IHiittemeister et aljt 1993a) are selected as intensity peaks in the 
CS maps from Ballv et al. ( 1987), therefore these correspond to 
high density regions, as CS is a dense tracer. Our sources, on 
the contrary, correspond to shock positions, therefore expecting 
that the amoun t of warm molecular gas w ill be larger than at the 
positions from lHuttemeister et alJ (Il993al) . 



4.2. On the heating and cooling of the molecular gas in the 
GMLs and in the CMZ 

Although the number of positions observed in the halo is rather 
limited, a remarkable result obtained in this work, is the single 
high kinetic temperature regime (> 95 K) for the halo sources, 
in contrast with the two temperature regimens (cool and warm) 
present throughout the CMZ. 

Therefore, the question that naturally arises is whether the 
cooling in the CMZ is more effective than in the halo positions, 
or if the heating mechanism in the GMLs is so efficient that the 
molecular gas has no time to be cooled down. 



4.2.1 . The cooling of the gas in the GC region 

In the following, we describe the cooling rates for H2 and CO 
emis sions, and for the ga s -dust coupling . 

iGoldsmith & Langerl (1 1978b derived the temperature depen- 
dence of the total cooling rate for a variety of molecular hy- 
drogen densities and at velocity gradient of 1 km s" 1 pc _1 for 
the temperature range of 10 to 60 K. According to them, for 
n(H2) < 10 3 ctrT 3 the cooling is dominated by CO; and for 

10 3 < n(H2) < 10 5 cm 3 , CI, O2, and the isotopic species of 
CO, contribute with the 30% to 70% of the total cooling. We 
use the expressions of the total cooling rates (A tota i) derived 
by them for the different physical conditions of the molecular 
gas i n the positions where this formulation is valid (see Table 
IA21 . iNeufeld et all (fl995h and lNeufeld & Kaufman! Jl99l de- 
rived the radiative cooling rates in a wider range of temper- 
atures (10 < T < 2500 K), and stated that the dominant 
coolants for the molecular interstellar gas are CO, H2, O, and 
H2O. The efficiency of the cooling depends of the gas temper- 
ature, the H2 particle density, and the optical depth parameter 
N(H2) = n(H2)/\dv/dz\, where dv/dz is the velocity gradient 
along the line of sight. The latter depends on the geometry and 
veloc ity structure assumed (see Table 1 in Neufeld & Kaufman 
119931) . Here we assume A>(H 2 ) = N(H 2 )/Av; where N(U 2 ) is 
the H2 column density estimated in column 1 1 of Table [4] and 
Av is the measured line width of CS(2-1) (Table POT l. We com- 
pute the total cooling rates according to their model in Table lA.2l 
by spline interpolation of the values for their molecular cooling 
function (Neufe ld etaD 19951). 

Later. ILe Bourlot et alJ (1 19991) derived the cooling rate by H2 
(Ah,) for a wider range of physical parameters (100 < T^ n < 

10 4 K; 1 < «(H 2 ) < 10 8 cm- 3 ). We used the FORTRAN subrou- 
tine provided by them only towards the sources allowed by the 
parameters range, i.e., for the warm gas (Table IA~2l . 

The cooling rate for the cou pling of the dust and gas is given 
bv fcoldsmith & Lange rll 19781) 



A 



gd 



2.4 x 10- 33 r g 1/2 



T l J £ (T g -T d )n 2 (H 2 )ergs- l cm- 3 



(1) 



where the grains parameters ( grain s i ze and a ccommodation co- 
efficien t) where taken from iLeungl (Il975l) . IHiittemeister et al] 
(fl993l) argued that the cold gas component is coupled to the 
dust temperature at high densities, therefore, the dust in the GC 
region would be a cooling agent. They proposed that the density 
of the cold gas must be at least an order of magnitude higher 
than that of the hot gas; otherwise the hot gas would also cooled 
down. We find that the gas density of the hot component is only 
slightly lower than in the cold component (a factor ~ 2 - 4 from 
the CS data in Table[3]l. Although t he gas-dust coup ling becomes 
significant at n{H2) ~ 10 5 cm" 3 (lJuvela & Ysardll201 ll) . which 
is reached only in the cool component regime in few positions 
in the CMZ (Table [3}, we estimated this cooling rate for all the 
positions (Table lA2l . 

Estimating tota l coo ling rates, e.g., following 
IGoldsmith & Langerl ([1978), one should account for two 
factors: the depletion of coolant species, and the lack of 
processed gas in the halo. The depletion of the coolant species 
can increase the gas temperature at low and moderat e densities 
of n(H 2) ^ 10 4 crrT 3 . This effect was studied by iGoldsmithl 
(2001) in dark clouds. In the physical conditions of the GC, 
the depletion of the coolant species is unlikely due to the high 
densities and low temperatures which are needed to have this 
effect. Second, as we noted before, in the density range of the 
GC clouds, CI, O2, and t he isotopes of CO, con tribute 30% to 
70% of the total cooling. Riauel me et al.l (1201 Oal) found a high 
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12 C/ 13 C isotopic ratio towards the disk-halo connection regions 
(halo, disk XI) which they interpreted as non-processed gas 
being accreted towards the GC. If this interpretation is correct 
(since the gas in the loop is less processed than the gas in the 
CMZ), one would expect a lower metallicity and, hence, a less 
efficient cooling by molecular or atomic lines. On the other 
hand, a high isotopic ratio was also found towards the disk XI 
positions, and those positions have indeed a cool temperature 
regime. 

The low temperature regime toward the XI orbit positions 
can be explained by t he large H2 cooling rate derived from 
iLe Bourlot et al.1 JT9 99) and the large total cooling rates derived 
from lNeufeld et al] (1 19951) (see Table lA.2"T >. However, the cooling 
rate in the Halo positions are also large, and those positions do 
not present the cool temperature regime. 

In the following, we estimate the heating rates affecting the 
GC region to study if a rise of the heating mechanism should be 
the responsible of the lack of low temperature regime in the halo 
sources. 

4.2.2. The heating of the gas in the GC 

We find high gas kinetic temperatures for basically all the ob- 
served positions of our sample. Therefore the heating mecha- 
nism responsible for the widespread high temperatures should 
apply for the gas in the entire GC region, with little effect on the 
dust. The heating mechanism in the GC should be different from 
that the heating of warm clouds (Tk; n > 100 K) in the disk, where 
the molecules are heated by collisions with hot dust, heated by 
embedded stars. 

Several heating mechanisms have been proposed to explain 
the high kinetic temperatures in the GC region. 

1. Cosmic rays heati ng: Heating by a la r ge flux of lo w en 



ergy cosmic rays (Giisten et al.l 1 198 it Morris et al. Il983t 
iHuttemeister et al l Il993at lYusef-Zadeh et all l2007al) . This 
mechanism requires a cosmic ray ionization rate ((cr) of one 
or two orders of magnitude larger than the Galactic value of 
10 _17 s _1 , which will influence also the gas-phase chemistry, 
increasing the atomic hydrogen due to the increased cosmic 
ray dissocia tion rate of H?, and also mo lecular ion emission 
like HCO + . lYusef-Zadehetail d2007allbh argue that such a 
high ionization rate is found in the GC region. From absorp- 
tion lines of H3 origin ating from a diffu se, hot molecular 
component in the CMZ. lGoto et alj d2008l) estimated an ion- 
ization rate of ~ 10~ 15 We assume a similar value in 
the gas observed by us, although it is presumably denser and 
cooler. The heating rate depends on n and £cr- It is difficult 
to derive it for each observed position, and this mechanism 
cannot be ruled out. 
2. X-ray heating: Heati ng by an extended diffuse source of 
soft X-ray emission (Watso n et al.l ll98"Th is unlikely since 
the X-ray emission is less extended than the NH3 emis- 
sion, and there is no obvious source for the required lu- 
mino sity in the extended soft X-ray emission (Morr is"et alJ 
19831). Heating by an ex tremely hot plasma emitting X-rays. 
Nagavama et alj d2007l) found that the NH3 emitting region 
in the CMZ, with 7^ = 20 - 80 K and T kin > 80 K, 
is surrounded, i n the longitude-ve locity space, by a high- 
pressure region (Sawadaetal. 2001), where the gas is less 
dense and hotter {n(H2) < 10 3 cm" 3 , T^ n > 100 K). Because 
the high pressure region is found to be coincident with the 
hot emitting X-rays, they argued that the thermal energy ra- 
diated from the hot plasma emitting X-ray plasma can heat 



the gas in the high pressure region. The heating rate due to 
this mechanism is very uncertain, but Ao et al., (2012, in 
prep.) argue that this mechanisms is not able to account for 
the high kinetic temperatures in the GC region. 

3. Ion-Slip heating: T he GC is pervaded b y a magnetic field 
of few mG (see e.g., iFerriere et al 112001 ." and their presence 
can also influence the heating of the GC. The ion-slip heating 
has been proposed for the molecular clouds in the GC region, 
where the heating rate depends on the magnetic field B, ion- 
ization rate, neutral number density n n , ion number density 
flu number of collision per second, the reduced mass of the 
ions and ne utrals, and the scale of the cloud R (in pc) (see, 
IScaldfl977l) . Because of the uncertainty of some on these 
values for each observed position, it is difficult to estimate 
the heating rate. 

4. Ultraviolet heating: The high NH3 abundances in the GC re- 
gion require effective shielding from UV radiation because 
ammonia is easily photodisociate d by ultraviolet radiation 
(Rodriguez-Fernand ez et al.1l200ll) . This is also confirmated 
by the large abundance of the HNCO molecule, which is also 
photodissociated by UV radiation. We discard UV heating in 
the GC region. 

HNCO could be formed via gas phase reactions, but for- 
mation in grains seems to be more efficient (see e.g., 
iMartfn et alJl2008l and references therein). Presumably this 
molecule is releas ed to the ISM by grain erosion and/or dis- 
ruption by shocks dZinchenko et al.ll2.000t) . and is easily pho- 
todissociated by UV radiation. The shocks that release the 
molecule from the grain mantles should be slow enough to 
not dissociated the molecule. We find that the X(HNCO) is 
low in the Disk XI -1 source, where we expect shocks, but no 
UV emission. Also we can see an enhancement toward the 
Halo 1 source (cloud number 2), and the DiskX2-l source, 
which could be due to the shocks present in these regions, are 
strong enough to evaporate the molecule from the grain man- 
tles but too slow to not dissociate the HNCO. This molecule 
can be used to trace the shocks properties throughout the GC 
region. 

5. Shocks: The dissipation of mechanical turbulence through 
shocks would offer a compelling answer, because the 
GC region shows an ubiquit ous presence of shocks, as 
traced by the SiO em ission (Martin- Pintado et al ] I2000L 
Riaue ime et al.ll2010bl) . The heating rate due to the dissipa- 
tion of tur bulence (of velocity v, on a spatial scale of R c ) is 
given by dBlackll 19871) 



3.5 x 10 



lpc 

~k7 



ergs 



(2) 



We derive the heating rate of the dissipation of turbu- 
lence for each position in Table IA.2I We use R c = 5 pc 
dRodrfguez-Fern andez et al. 1 12001 ). It is important to note 
that this equation is highly dependent of v t , therefore if the 
molecular clouds are not resolved in the beam size, the heat- 
ing rate would be overestimated. 

The origin of this turbulence in the positions studied in this 
work, can be the follow ing: the large sca le dynamics in a 
barred potential model dBinnev et al.ll 199 II) can produce the 
shocks found towards the high velocity clouds associated 
with the XI orbits (in the 1.3 complex, and SgrC). This is 
supported by the higher kinetic temperatures found in the 
diskXl sources than in the diskX2 sources (Table©. In our 
"halo" positions the shocks can be produced by the GMLs 
scenario, which is supported by the broad velocity width at 
the foot point of the loops. However, high kinetic temper- 
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atures and the large SiO emission are widespread through- 
out all the GC region. Supernova or hypernova explosions 
could cause for the high temperature found in the lower 
velocity c omponents (in our notation, diskX2) of the 1°.3 
complex (Tana ka et al.l 12007b . For the lower longitudes in 
the CMZ, the heating could be explained by interactions 
with SNRs close to SgrA and also with non-thermal fila- 
ments in the radio arc, and cloud-clou d collisions, and ex- 
pandin g bubbles in the vicinity of Sgr B Mart in-Pintado et all 
( 1997). Cl oud-cloud collis ions were also proposed in the GC 
region dWilson et al.lll982l) . which is favored by the la rge 
linewidfh typical in the GC clouds. iGiisten et al.l (1 19851) ar- 
gue that if this mechanism is acting in the molecular gas, 
the linewidfh and t he te mperatures should be correlated. 
Mauersbe rger et al.l (Il986al) found such a correlation for the 
clouds obse rved in the GC, which woul d support this mech- 
anism. Like [Hiittemeiste r et al.l (fl99 3a). we do not confirm 
such a clear correlation between the rotational temperature 
associated with the inversion transitions (2,2), (4,4), and 
(5,5) (Table |2]i and the Doppler width of our sources (see 
Fig 0. For the few observed positions we can neither con- 
firm nor reject the cloud-cloud collisions as primary heating 
mechanism in our sources. 

The heating rates derived from Eq. [2] for the halo and disk po- 
sitions are similar. Therefore, it is probably that this mecha- 
nism by itself is not causing the lack of the low temperature 
regimen observed in the halo positions. An extra heating input 
would be required for the halo positions. iTorii et ai1 d2010a) pro- 
pose that the gas in the foot point of the GML is heating by 
C-shocks, and that the warmest region of the foot point (the 
"U shape") would also be heated by magnetic reconnection or 
by upward flowing gas bounced by the narrow neck in the foot 
point. They estimated that the total available energy (considering 
the magnetic and gravitational energy) injected to the U shape is 
1.8 - 2.6 x 10 37 ergs~' in ~ 10 6 years . Consid er a cloud size of 3 
pc radius, as estimated bv lTorii et al.l (l2010al) . the heating rate is 
T ~ 3.9 -7.8 x 10~ 21 erg crrT 3 s . This heating rate is negligible 
in comparison with the values obtained from the dissipation of 
turbulence alone. Alternatively, the CMZ can be understood as 
a highly turbulent medium, where many phenomena are taking 
place (shocks produced by Galactic potential, SN explosions, 
star formation, interaction with non-thermal filaments, cosmic 
rays, the presence of a supermassive black hole, etc) and coex- 
ist in the central few hundred parsec of the Galaxy. All of these 
phenomena modify the physical parameter of the region. The 
two kinetic temperature regimes present in this reg ion, which are 
in pressure equilibrium (Htitteme ister et al.ll 19981) . could be the 
results of the interplay of the different phenomena mentioned 
above. On contrary, in the GMLs scenario, the gas goes down 
towards the Galactic plane following the magnetic field lines, 
which lead a tidy movement of the gas producing shocks front at 
the foot point of the loops. The continuous shock at the foot point 
is not affected by other phenomena (like, e.g., star formation), 
therefore the gas is continuously heated. High resolution maps 
of the foot point of the GMLs are needed not only to confirm 
this hypothesis but also to resolve the linewidth of the molecular 
clouds to better estimate the heating rate for dissipation of tur- 
bulence (Eq. |2j. Summarizing, in spite of the limited number of 
positions studied in this work, we propose that the high kinetic 
temperature found in all of the sources are produced by shocks. 
We discarded UV heating due to the large abundance of HNCO 
and NH 3 molecule. From our data it is not possible to confirm or 
rule out x-ray and ion-slip heating. 
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Fig. 3. T ro , is plotted versus the average velocity linewidth from 
all the metastable inversion transitions. The T rot correspond to 
the values in bold from the Table [2] 



5. The ammonia ortho-to-para ratio and its 
implications on kinetic temperatures now and 
there 

Radiative and collisional transitions between ortho-NH3 (K = 
0,3,6,9...) and para-NH 3 (K = 1,2,4,5...) are forbidden be- 
cause of their different nuclear spins. The time scale of conver- 
sion between ortho a nd para species is 10 6 yrs in the gas phase 
dCheung et al.lll969l) . Therefore orfho-NH 3 and para-NH 3 can 
be almost treated as different species. The spin temperatures be- 
tween ortho and para-NH 3 , or the ortho/para ratio, may reflect 
the conditions at the time of the formation of NH 3 , while rota- 
tional temp_eraturewifli^tte present condi- 
tions dHo & Towneslll983HUmemoto et al.lll999l) . 

A ratio of 1 .0 would be expected if NH 3 is formed in the gas 
phase reactions at high temperatures, and a larger value could be 
explained by a formation on cold dust grains and a subsequent 
release into the gas phase. For temperatures of about 10 K, the 
typical temperature of cold dust, the ortho-to-para ratio will be 
larger than 2. 

iDulieul d201 ll) . however, points out that there might be a 
differential desorption on ortho and para molecules from dust 
grains and that this, and not the formation process, determines 
the gas phase ortho/para ratios. 

From the LVG analysis of the ammonia molecule, we de- 
rive an ortho-to-para ratio average for all of our sources of 0.70 
with a standard deviation of 0.20 (see Table |7). In this cal- 
culation, we take into account only the warm component of 
the para-NH 3 with the ortho-NH 3 . While the ortho-NH 3 tran- 
sitions are probably tracing higher kinetic temperatures than the 
para-species (see Table IA. 31 . the ortho-to-para ratio derived us- 
ing RADEX for the warm kinetic temperature component ap- 
proaches to the statistical equil i brium v alue. This result is con- 
sistent with iHiittemeister et al.l (|T993a) estimation of the (3,3) 
column density. Non thermal emission has been pr edicted for the 
NH 3 (3,3) line by Walmslev & Ungerechts ( 1983), and has been 
observed e.g. by [Mauersberg er et al.l (Il986a | b|). Also the (6,6 ) 
line might be a maser source in some case ( Lebron et al.ll201 lb . 
Our RADEX calculations were made assuming a background of 
2.7 K. We have verified that in our models the ortho lines tend 
to be weak masers for a large part of our parameter space. 
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Table 7. Ortho-to-para NH 3 ratio from LVG model 



Source 




AT/ \TIT \I\T/ A T T T \b 

N(o - NH 3 )/N(p - NH 3 ) 


Halol 


1 


0.83 




2 


0.63 




i 


0.05 


Halo4 


1 


0.48 


DiskXl-1 


1 


1.0 


DiskX2-l 


1 


0.76 


DiskXl-2 


1 


0.93 


DiskX2-2 


1 


0.76 


Diskl 


1 


2.02 c 




2 


0.72 c 


Disk2 


1 


0.59 




2 


0.63 




3 


0.34 



" Cloud number defined by the different velocity components. 

" We take into account only the warm component of the para-NH 3 

with the ortho-NH 3 (see text for details). 

c This value was not considered in the average of the ortho-to-para 
ratio, because it is possible that this position does not have the high 
kinetic temperature regimen. 



6. Conclusions 

1. We have used the metastable inversion transitions of NH3 
from (J, K) = (1,1) to (6,6) to derive the gas kinetic tem- 
perature toward six positions selected in the Galactic cen- 
tral disk and for higher latitude molecular gas. We also ob- 
served other molecules like SiO, HNCO, CS, C 34 S, C ls O, 
and 13 CO, to derive the densities and to trace different physi- 
cal processes (shocks, photodissociation, dense gas) which 
was used to reveal the heating mechanisms affecting the 
molecular gas in these regions. 

2. The GC molecular gas consists of roughly two kinetic tem- 
perature components in the CMZ. Only the warm kinetic 
temperature regime is found in the "halo" positions, and in 
the Disk 2 position, which corresponds to SgrB2. The results 
obtained in this paper apply for the disk and halo positions 
defined in this work within the GC region, and do not repre- 
sent the conditions for the Galaxy as a whole. 

3. The kinetic temperatures are high, not only in the typical GC 
clouds, but also in the high latitude and high velocity clouds 
observed in this paper. 

4. Shocks are a compelling heating mechanism of the molec- 
ular clouds in the GC region. This is supported by the high 
gas kinetic temperature and by the increased SiO abundance 
in the location where shock are expected. Due to the fragile 
nature of the HNCO (enhanced by shock but easily photodis- 
sociated by UV radiation and strong shocks), this molecule 
could be used to reveal the characteristic of the shock. Other 
heating mechanisms previously proposed for the GC, how- 
ever, cannot been ruled out. 

5. The high kinetic temperatures found in both, the XI orbits 
and in the foot point of the GMLs, seem to support to the 
large scale dynamics induced by the bar potential, and the 
GMLs scenario as origins for the shocks respectively. 
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Table A.l. Fractional abundances of NH 3 with respect to SiO, HNCO, CS, C 34 S, and of SiO and HNCO with respect to CS and C 34 S 



Source 


rf 




ST7T7T7 — \ 






N(ii lU) 


A>(i IV) 


N(HNCO) 


N(HNCU) 


N(ii lU) 




N(SiO) 


N(HNCO) 


N(CS) 


N(C M S) 


N(CS) 


N(C 34 S) 


N(CS) 


N(C 34 S) 


N(HNCO) 


Halol 


1 


11.97 ± 3.64 


53 64 + 16 28 


3.99 ± 1.21 




> 0.33 




> 0.07 




> 4.48 




2 


9.11 +2.68 


14.59 ± 3.09 


2.32 ± 0.41 


33.40 ± 7.84 


0.25 ± 0.09 


3.67 ± 1.35 


0.16 ±0.04 


2.29 ± 0.70 


1.60 ± 0.57 




3 


3.41 ±2.67 




1.01 ±0.76 




0.29 ±0.31 










Halo2 


1 

2 










0.14 ±0.15 
> 0.12 


1.69 ±2.01 


0.29 ±0.31 


3.47 ± 4.05 


0.49 ± 0.53 


Halo3 


1 

2 










0.17 ±0.09 


2.60 ± 1.50 


0.20 ±0.10 
0.33 ±0.17 


3.08 ± 1.58 
3.89 ±2.58 


0.85 ± 0.42 


Halo4 


1 


32.32 ± 18.12 


8.60 ± 3.46 


5.46 ± 2.84 


81.43 ±49.50 


0.17 ±0.12 


2.52 ± 1.97 


0.64 ± 0.38 


9.47 ± 6.41 


0.27 ± 0.17 


Halo5 


1 










0.19 ±0.07 


2.18 ±0.80 


0.49 ±0.15 


5.60 ± 1.64 


0.39 ±0.14 


DiskXl-1 


1 


8.91 ± 2.30 


12.46 ± 2.51 


2.69 ± 0.65 


52 80 + 13 87 


30 + 08 


5.93 ± 1.74 


0.22 ± 0.05 


4.24 ± 1.04 


1.40 ± 0.34 




I 




5.57 + 1.70 


2.48 + 1.04 


SR 46 + ?4 9? 


17 + 010 


4 08 + ? 44 


44 + ?? 




3Q + o ?0 


DiskXl-2 


1 


14.32 ± 2.76 


9.83 ± 1.76 


2.78 ± 0.48 


42.41 ± 6.94 


0.19 ±0.04 


2.96 ± 0.62 


0.28 ± 0.06 


4.31 ±0.86 


0.69 ± 0.15 


DiskX2-2 


1 


111.3 ±34.7 


14.52 ±3.28 


2.60 ± 0.70 


29.67 ± 7.91 


0.02 ± 0.01 


0.27 ± 0.08 


0.18 ±0.04 


2.04 ± 0.47 


0.13 ±0.04 


Diskl 


1 


6.67 ± 0.00 


> 2.84 


> 1.11 


93.33 ± 0.00 


> 0.17 


> 14.00 


> 0.39 


> 32.83 


> 0.43 




2 


17.28 ±6.91 


4.59 ± 1.32 


1.95 ±0.72 


45.96 ± 18.58 


0.11 ±0.06 


2.66 ± 1.51 


0.42 ± 0.20 


10.02 ± 4.98 


0.27 ±0.13 


Disk2 


1 


74.31 ±23.50 




24.65 ± 7.26 


71.03 ±20.94 


0.33 ±0.14 


0.96 ± 0.41 










2 


89.23 ±0.19 


> 1.90 


30.43 ± 0.06 


88.60 ±0.19 


>0.34 


>0.99 


> 16.01 


> 46.63 


> 0.02 




3 






14.30 ± 1.86 















" Cloud number defined by the different velocity components. 
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Table A.2. Estimation of heating and cooling rates for each source 



Source 


n" 


Tkin 

[K] 


n(H 2 ) 
[xl0 4 cnr 3 ] 


A be 
A 

total 


\otal 


A 


A . 


^ lurb 


Halol 


1 


115 


1.00 




14.0 


0.78 


0.018 


27.5 


Halol 


2 


90 


8.49 




242 


4.36 


0.86 


58.4 


Halol 


3 


135 


1.58 




46.8 


31.4 


0.057 


12.9 


Halo2 


1 


113 


1.26 




21.6 


0.015 


0.028 


31 


Halo2 


2 


113 


0.97 




13.3 


0.82 


0.016 


16.5 


Halo3 


1 


113 


1.50 




30.7 


1.3 


0.039 


62.6 


Halo3 


2 


113 


3.55 




136 


3.1 


0.22 


178 


Halo4 


1 


95 


2.75 




44.9 


1.4 


0.10 


41 


Halo5 


1 


95 


7.47 




210 


3.9 


0.74 


51.4 


Disk-Xl 


1 


38 


13.0 


1.5 


36.7 




0.081 


849 






300 


3.55 




534 


261 


0.79 


232 


Disk-X2 


1 


38 


5.62 


1.5 


12.7 




0.015 


234 






100 


2.70 




48.5 


1.2 


0.11 


113 


Disk-Xl 


2 


52 


13.3 


3.6 


101 




0.51 


76.3 






215 


5.62 




746 


112 


1.3 


32.2 


Disk-X2 


2 


28 


14.1 


0.60 


16.5 




0.57 


158 






95 


6.20 




150 


3.3 


0.51 


69.4 


Diskl 


1 


23 


1.12 


0.071 


1.26 




0.005 


5.29 






154 


0.16 




1.37 


0.51 


0.001 


0.76 


Diskl 


2 


38 


4.57 


0.28 


25.8 




0.010 


31.4 






82 


2.51 




45.1 




0.064 


17.3 


Disk2 


1 


68 


4.97 


1.3 


48.6 




0.17 


250 






200 


2.47 




179 


34 


0.23 


124 


Disk2 


2 


145 


2.51 




55.8 


8.1 


0.16 


32.5 


Disk2 


3 


50 


15.8 


3.2 


86.8 




0.60 


129 






80 


11.6 




224 




1.3 


94.7 



8 Cloud number denned by the different velocity components. 



* from lGoldsmifh & LaneeJ jl978h. 

c from lNeufeld et alj jl995h ; Neufeld & Kaufmarj d!993h 
d using program provided by Le Bourlot et al.Ul999h . 

f [xlO^'ergs-'cnr 3 ] 

Table A.3. Column density derived for NH3 (3,3) and NH3 (6,6) using RADEX 



Source 


a 


n(H2) 


T(kin) 


N(ortho-NH3) 


Halol 


1 


1.00 x 10 4 


> 300 


1 x 10 14 


Halol 


2 


8.49 x 10 4 


> 300 


1.26 x 10 14 


Halol 


3 


1.58 x 10 4 


> 300 


0.45 x 10 14 


Halo4 1 


2.75 x 10 4 


> 170 


0.95 x 10 14 


DiskXl-1 1 


3.55 x 10 4 


> 300 


1.58 x 10 14 


DiskX2-l 1 


2.70 x 10 4 


> 200 


2.51 x 10 14 


DiskXl-2 1 


5.62 x 10 4 


> 280 


1.41 x 10 14 


DiskX2-2 1 


6.20 x 10 4 


265 


0.76 x 10 14 


Diskl 




0.16 x 10 4 


> 300 


0.32 x 10 14 


Diskl 


2 


2.51 x 10 4 


> 300 


0.32 x 10 14 


Disk2 


1 


2.47 x 10 4 


200 


8.32 x 10 14 


Disk2 


2 


2.51 x 10 4 


> 300 


12.6 x 10 14 


Disk2 


3 


11.6 x 10 4 


130 


4.79 x 10 14 



Cloud number denned by the different velocity components. 
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Table A.4. Results from Gaussian fit and optical depth from NH 3 method. 



Source 


Cloud 


Transition 


VLSR 


Tmb 


Avi/2 


Integrated intensity 15 


rms 


T 




number 


(j,k) 


r 1 -1 1 

[km s '] 


[K] 


r i -In 

[km s '] 


[K km s~'] 


[mk] 




Halol 


1 


(1,1) 


85.6 ±0.4 


0.10 


29.0 ± 0.4 


3.19 ±0.04 


52 


< 0.40 




2 


(1,1) 


115.8 ±0.4 


0.32 


23.6 ± 0.4 


7.91 ±0.04 


52 


0.10 ±0.09 




3 


(1,1) 


139.3 ±0.4 


0.09 


18.8 ±0.4 


1.77 ± 0.04 


52 


< 0.64 




1 


(2,2) 


83.1 ± 0.4 


0.07 


37.8 ±0.4 


2.94 ± 0.04 


51 


< 0.86 




2 


(2,2) 


116.7 ±0.4 


0.35 


18.4 ±0.4 


6.77 ± 0.04 


51 


0.21 ±0.01 




3 


(2,2) 


147.1 ±0.4 


0.06 


35.5 ± 0.4 


2.29 ± 0.04 


51 


4.9 ± 2.0 d 




1 


(3,3) 


87.7 ± 0.4 


0.22 


38.8 ±0.4 


9.23 ± 0.07 


51 


0.10 ±0.03 




2 


(3,3) 


117.2 ±0.4 


0.85 


16.8 ±0.4 


15.16 ±0.07 


51 


0.10 ±0.03 




3 


(3,3) 


138.6 ± 0.4 


0.19 


20.8 ± 0.4 


4.29 ± 0.07 


51 


0.16 ± 0.01 




1 


(4,4) 


82.0 ± 0.4 


0.06 


32.9 ± 0.4 


2.03 ± 0.02 


50 






l 


(4,4) 


1 1 O H i A A 

lis./ ± 0.4 


0.23 


H 1 C | A A 

21. D ± U.4 


5.34 ± 0.02 


CA 

50 






3 


(4,4) 


137 .4 ± 0.4 


0.05 


16.6 ± 0.4 


0.91 ± 0.02 


50 






1 


(5,5) 


87.7" 


0.12* 


38.8" 


0.467* 


41 






2 


(5,5) 


115.7 ± 0.4 


0.15 


12.4 ± 0.9 


1.92 ± 0.13 


41 






3 


(5,5) 


131.0 ± 0.9 


0.06 


11.9 ± 3.2 


0.74 ± 0.09 


41 






1 


(6,6) 


87.7" 


0.13* 


38.8" 


0.613* 


41 






f\ 

z 


(6,6) 


1 i a c i n o 

120.5 ± 0.8 


0.30 


17.3 ±2.1 


C AO i t\ C\ 

5.48 ± 0.61 


41 






3 


(6,6) 


138.5 ± 1.9 


0.12 


14.4 ± 2.9 


1 TO f\ C A 

1.78 ± 0.54 


41 




Halo 4 


1 


(1,1) 


196.4 ±0.9 


0.22 


30.2 ± 2.1 


6.99 ± 0.42 


83 


< 0.31 




1 


(2,2) 


194.7 ± 1.1 


0.19 


34.1 ± 2.9 


6.74 ± 0.47 


84 


1.72 ± 0.53 




1 


(3,3) 


195.5 ± 0.5 


0.37 


20.3 ± 1.2 


7.93 ± 0.38 


80 


< 0.31 




1 


(4,4) 


194.5 ± 1.1 


0.09 


16.0 ± 2.9 


1.44 ± 0.20 


49 






1 


(5,5) 


195.5" 


0.15* 


20.3° 


0.417 


50 






1 


(6,6) 


199.4 ± 1.7 


0.09 


24.3 ± 4.0 


2.26 ± 0.32 


51 




DiskXl-1 


1 


(1,1) 


183.1 ± 1.0 


0.20 


48.1 ± 2.2 


10.43 ± 0.45 


74 


0.86 ± 0.27 




1 


(2,2) 


182.9 ± 1.3 


0.15 


32.1 ± 2.6 


5.10 ± 0.38 


71 


0.10 ± 0.11 




1 


(3,3) 


184.1 ±0.4 


0.43 


34.8 ± 1.0 


15.88 ±0.39 


71 


0.10 ±0.02 




1 


(4,4) 


188.1 ± 1.1 


0.14 


24.7 ± 2.2 


3.66 ± 0.33 


74 






1 


(5,5) 


187.1 ± 1.7 


0.09 


36.4 ±4.1 


3.56 ± 0.34 


63 






1 


(6,6) 


188.3 ± 1.1 


0.17 


35.0 ± 2.6 


6.25 ± 0.41 


59 




DiskX2-l 


1 


(1,1) 


96.4 ± 0.4 


0.48 


40.8 ± 1.0 


20.91 ± 0.42 


74 


< 0.84 




1 


(2,2) 


96.5 ± 0.5 


0.32 


30.7 ± 1.4 


10.36 ± 0.38 


71 


0.10 ±0.02 




1 


(3,3) 


98.1 ± 0.2 


0.74 


32.9 ± 0.6 


25.95 ± 0.39 


71 


0.10 ±0.01 




1 


(4,4) 


96.0 ± 0.6 


0.19 


15.7 ± 1.7 


3.16 ±0.28 


74 






1 


(5,5) 


AO e~ i 'i 'i 

98.6 ± 2.2 


a nn 

0.09 


A 1 C i A C 

41.5 ± 4.0 


o A 1 i a n 

3.91 ± 0.37 


63 






1 


(6,6) 


103.4 ± 1.1 


0.17 


30.8 ± 2.6 


5.70 ± 0.40 


59 




DiskXl-2 


1 


(1,1) 


69.1 ± 0.5 


0.32 


22.7 ± 1.3 


7.77 ± 0.37 


78 


0.10 ±0.03 




1 


(2,2) 


69.4 ± 0.4 


0.32 


14.5 ± 1.3 


4.99 ± 0.33 


80 


0.10 ±0.06 




1 


(3,3) 


69.1 ± 0.2 


0.87 


17.3 ± 0.6 


16.10 ± 0.42 


91 


0.10 ±0.04 




1 


(4,4) 


69.9 ± 0.4 


0.25 


12.2 ± 1.4 


3.23 ± 0.26 


71 






1 


/c c\ 
(j,j) 


£~H A i a H 
D/.U ± U. / 


U. 13 


1 "2 1 i 1 A 

13.1 ± 1.4 


1 O A i A OA 

1.84 ± U.2U 


/CO 

oz 






1 


(6,6) 


74.1 ± 0.6 


0.24 


17.8 ± 1.7 


4.58 ± 0.33 


61 




Disk X2-2 


1 


(1,1) 


-50.2 ± 0.8 


0.25 


37.9 ± 1.9 


10.22 ± 0.46 


78 


1.00 ± 0.29 




1 


(2,2) 


-53.0 ± 1.2 


0.15 


22.1 ± 3.1 


3.51 ± 0.39 


80 


< 0.28 




1 


(3,3) 


-47.1 ± 0.6 


0.35 


22.9 ± 1.4 


8.41 ± 0.45 


91 


< 0.23 




1 


(4,4) 


-54.1 ± 0.6 


0.14 


6.7 ± 1.5 


0.99 ±0.18 


71 






1 


(5,5) 


-47.1° 


0.19* 


22.9" 


0.55* 


62 






1 


(6,6) 


-49.3 ± 2.8 


0.07 


21.0 ±4.9 


1.31 ± 0.31 


61 




Disk 1 


1 


(1,1) 


64.5 ± 0.4 


0.08 


20.0 ± 0.4 


1.73 ± 0.03 


61 


< 1.25 




2 


(1,1) 


78.2 ± 0.4 


0.16 


18.4 ± 0.4 


3.05 ± 0.03 


61 


< 0.34 




1 


(2,2) 


59.3 ± 0.8 


0.08 


5.1 ± 1.6 


0.45 ±0.15 


63 


< 1.35 




f\ 

z 


(2,2) 


73.4 ± 0.5 


a 1 *c 

0.16 


8.5 ± 1.5 


1 A *7 i A OA 

1.4/ ± 0.20 


63 


< 0.26 




i 
i 


1"? X\ 


JO. J ± l.U 


n 1 1 

U. 1 1 


1 1 7 j. 1 1 

ii./ ± j. i 


1 T.fi -a- 97 
l.jO ± u.z / 


oo 


<. U.JO 




z 


(J,J) 


hh a _i_ n A 
/Z.4 ± U.O 


a h i 
U.zl 


1/1 O _i_ A ^ 

14. y ± U.o 


i a 1 _i_ n (Yi 
j AL ± 0.03 


00 


< U.JZ 




1 


1 A A\ 

(4,4) 


C£fi 

3D 


U.lO 


1 1 Ha 
11. 1 


A lib 

0.33 








1 


(4,4) 


72" 


a i cb 

0.16 


15.0" 


A Q ob 

0.38 


53 






1 


(5,5) 


56° 


0.14* 


11.7" 


0.29^ 


47 






z 


is. ca 


1 Z 


U. 14 


1 J.U 




AH 
4 / 






1 


(6,6) 


56" 


0.15* 


11.7" 


0.39* 


50 






2 


(6,6) 


72° 


0.15* 


15.0" 


0.44* 


50 




Disk 2 


1 


(1,1) 


46.1 ±0.1 


2.2 


24.1 ±0.3 


56.34 ± 0.65 


94 


0.10 ±0.10° 




2 


(1,1) 


73.9 ±0.1 


3.0 


18.4 ±0.2 


58.21 ±0.40 


94 


3.00±0.10 c 




3 


(1,1) 


90.2 ±0.1 


3.6 


16.7 ±0.1 


64.73 ± 0.33 


94 


0.10 ±0.10° 




1 


(2,2) 


44.9 ±0.1 


1.8 


22.6 ± 0.4 


42.52 ±0.58 


106 


0.10 ±0.01 




2 


(2,2) 


74.3 ± 0.2 


3.2 


17.2 ±0.3 


57.55 ± 1.43 


106 


0.15 ±0.01 
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Table A.4. continued. 



Source Cloud 


Transition 






AVi/2 


Integrated intensity 6 


rms 


T 


number^ 


(i,k) 


[km s" 1 ] 


[K] 


[km s" 1 ] 


[Kkms~'] 


[mk] 




3 


(2,2) 


O/i 1 i n i 
y4.1 ± U.3 


H 1 

2. i 


i o t~ i n c 
is.o ± V.J 


a i no , 1 Q 1 

41. zo ±1.31 


1UO 


n i n i n n 1 
U.1U ± U.U1 


t 
1 


(3,3) 


4o.o ± (J. 1 


1 A 

3.4 


HO "7 i A 1 

28.7 ± 0.3 


1 ai /t , a o 

103. o ± 0.8 


no 

98 


U.1U ± 0.10 


2 


(3,3) 


HA C i A 1 

/4.J ± U. i 


1 A 

i4 


i a a i n 1 
14.4 ± U. 1 


Tit 1 , AT 

210.2 ± 0. / 


no 


ft in , n 1 nc 
U. 11) ± U. 11) 


3 


(3,3) 


95.8 ±0.1 


2.7 


18.6 ± 0.2 


<i . A C 1 

02.42 ± U.M 


98 


n 1 n i n 1 nc 

(J. ID ± U.10 


1 


(4,4) 


Jl). 1 ± U.O 


n o 


il.j ± 1.3 


20.23 ± u.y / 


HA 




9 

z, 


(A A\ 


ia n + n i 

/H-.W ± U.l 


.3. 1 


i s n + n ? 


T-7.JJ It l.UJ 


74 




3 


(4,4) 


96.2 ± 0.6 


0.48 


20.0 ± 1.2 


10.24 ± 0.57 


74 




1 


(5,5) 


53.5 ±0.1 


0.40 


35.2 ±0.5 


14.90 ± 0.20 


57 




2 


(5,5) 


73.0 ±0.1 


2.0 


15.9 ±0.1 


33.44 ± 0.28 


57 




3 


(5,5) 


97.3 ±0.6 


0.18 


16.3 ±1.4 


3.10 ±0.24 


57 




1 


(6,6) 


56.5 ± 0.3 


0.69 


15.9 ±0.8 


11.80 ±0.53 


51 




2 


(6,6) 


75.7 ±0.1 


2.4 


15.2 ±0.2 


38.76 ± 0.48 


51 




3 


(6,6) 


96° 


0.15* 


18.6" 


0.51* 


51 





a no detected. Value taken as reference from the (3,3) transition. 
* upper limits. 

c method NH3 did not converge properly. 

d the fit is not good, due to the no convergence of the method. 

' f T M Bdv. 

f In each source, different velocity components define different clouds. 
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Table A.5. Gaussian fits and column densities for SiO, C 34 S, HNCO 



Source 


Species 


Velocity Center 


A v 


T A 


J T K dV 


rms 


N 






LSR [km s- 1 ] 


[km s" 1 ] 


[K] 


[K kms- 1 ] 


[mK] 


[cm" 2 ] 


Halo 1 


SiO (2-W 


88 + 1 


33.5 ± 1.0 


0.096 


3.41 ± 0.01 


4.5 


S 73 + 01 x 10 12 






117.8 ± 0.1 


17.9 ±0.1 


0.993 


i 8 of, + o 08 


4.5 


31 82 + 14 x 10' 2 






1 33 5 + o 3 


20.7 + 0.7 


0.174 


J.O J IE \J.\Ja. 


4.5 


(S 4.^ + 0^ x 10 12 




olw(y 1 ) 


110 7+06° 

117./ IE U.U 


17 7+1 3 

1 / .Z. IE I. J 


U.UO J 


1 + 1 1 

1 .J J It U. 1 1 


f\ ^ 

\j.z> 


? (\f\ + 1 Q x 1 12 

A.\J\J IC U. 17 A 1U 




30 SiO (2-Y) 


1 20 + 8" 

J.Z.V.V -!- w • (J 


17.4 + 1.9 


0.063 


1.16 ± 0.11 


5.9 


2 04 + 19 X 10 12 




CS (2-1) 


84.8 ± 0.5 


27.5 ± 1.2 


0.211 


6.2 ± 1.8 


0.2 








116.5 ±0.1 


17.1 ±0.1 


1.627 


29.7 ± 0.2 


0.2 










18 2 + 06 

1 0< _!_ V/ ■ vj 


0.364 


7 06 + 09 

/ .\J\J _!_ \J.\Jjr 


0.2 






CS (3-2) 


83.1 ±4.1 


23.3 ±4.1 


0.075 


1.86 ± 0.96 


9.6 








116.5 ±4.1 


15.4 ±4.1 


1.517 


24.80 ± 0.96 


9.6 








130.7 ±4.1 


11.2 ±4.1 


0.279 


3.34 ± 0.96 


9.6 






C 34 S (2-1) 


















117. ± 0.3 


18.8 ± 0.8 


0.138 


2.75 ± 0.10 


5.6 


9.51 ± 0.34 x 10 12 




HNCO 


83.2 ± 2.0 


15.2 ±4.0 


0.023 


0.37 ± 0.09 


5.2 


10.8 ± 2.6 x 10 12 






117.2 ±0.4 


21.7 ± 1.1 


0.119 


2.75 ±0.11 


5.2 


81.0 ± 3.4 x 10 12 




C 18 O(l-0) 


89.0 ± 0.9 


26.8 ± 1.9 


0.057 


1.62 ±0.10 


4.7 


1.11 ±0.07 x 10 15 






115.9 ±0.2 


18.8 ±0.6 


0.108 


2.17 ±0.04 


4.7 


1.49 ± 0.03 x 10 15 






140.7 ± 1.4 


16.6 ±3.8 


0.022 


0.38 ± 0.07 


4.7 


0.26 ± 0.05 x 10 15 




13 CO(1-0) 


85.0 ± 0.4 


37.6 ± 1.4 


0.344 


13.77 ±0.4 


12.0 


0.94 ± 0.03 x 10 15 






114.2 ±0.2 


13.3 ± 0.4 


0.628 


8.92 ± 0.5 


12.0 


0.61 ± 0.03 x 10 15 






133.7 ± 1.4 


27.2 ±3.1 


0.167 


4.84 ± 0.6 


12.0 


0.33 ± 0.04 x 10 15 






169.5 ±2.3 


23.6 ± 6.0 


0.049 


1.23 ± 0.26 


12.0 


0.08 ± 0.02 x 10 15 


Halo 2 


SiO (2-11 

v. ? ivy / 


-84.7 ± 1.1 


25.5 ± 0.8 


0.046 


1.25 ± 0.01 


4.0 


2.10 ± 0.01 x 10 12 






-58.5 ± 1.0 


22 6 + 2 8 


0.045 


1.07 ± 0.11 


4.0 


1 80 + 19 x 10 12 




29 SiOf2-l^l 


-85* 


25* 




< 0.26 


6.8 


< 0.453 x 10 12 






-58* 


23* 




< 0.26 


6.8 


< 453 x 10' 2 




30 SiO (2-1) 


-58.6 ± 2.1" 


18.9 ±4.8 


0.020 


0.40 ± 0.07 


3.4 


0.71 ± 0.13 x 10 12 




CS (2-1) 


-80.2 ± 6.1 


25.8 ± 6.1 


0.188 


5.15 ± 0.30 


6.2 








-51.9 ± 6.1 


22.9 ±6.1 


0.192 


4.69 ± 0.30 


6.2 






CS (3-2) 


-77.6 ± 1.2 


22.2 ± 2.7 


0.091 


2.16 ±0.30 


6.5 








-51.8 ± 1.5 


22.3 ± 2.5 


0.078 


1.85 ±0.25 


6.5 






C 34 S (2-1) 


-81.7 ± 13.5 


47.8 ± 23.2 


0.006 


0.47 ± 0.24 


3.8 


1.63 ± 0.84 x 10 12 






-35.0 ± 14.9 


35.6 ± 21.9 


0.009 


0.24 ± 0.25 


3.8 


0.83 ± 0.85 x 10 12 




HNCO 


-70.4 ± 3.9° 


30 6 + 8 2 


0.032 


1 05 + 26 


5.1 


31 + 7 7 x 10 12 


Halo 3 


SiO (2-W 


-62.9 ± 1.2 


35.9 ± 3.4 


0.194 


7.43 ± 0.52 


8.4 


12 5 + 9 x 10 12 

1—.,..' _L_ \J * J /\ 1 V / 




29 SiOC2-11 


-56.6 ± 3.7 


23.1 ± 8.6 


0.015 


0.37 ±0.11 


4.8 


64 + 18 x 10 12 




30 SiO (2-1 "1 


-51.0 ± 8.2 


34.3 ± 15.5 


0.020 


73 + o 31 


3.3 


1 29 + 55 x 10 12 




CS (2-1) 


-64.7 ± 0.3 


31.4 ±0.8 


0.670 


22.42 ± 0.46 


8.2 








-13.9 ± 1.1 


33.0 ± 2.6 


0.186 


6.52 ± 0.44 


8.2 






CS (3-2) 


-63.6 ± 0.3 


18.6 ± 1.1 


0.513 


10.13 ± 0.57 


14.4 








-21.6 ±4.5 


42.7 ± 9.5 


0.077 


3.50 ± 0.60 


14.4 






C 34 S (2-1) 


-63.1 ± 1.2 


28.8 ±3.6 


0.051 


1.58 ±0.15 


6.5 


5.46 ± 0.52 x 10 12 






-12.6 ±2.4 


21.4 ±4.8 


0.023 


0.52 ±0.11 


6.5 


1.81 ±0.39x 10 12 




HNCO 


-67.3 ± 0.8 


35.3 ± 2.0 


0.083 


3.10 ± 0.15 


6.0 


91.5 ±4.3 x 10 12 






-10.6 ± 0.9 


21.7 ± 2.2 


0.055 


1.28 ± 0.11 


6.0 


37.7 ± 3.3 x 10 12 


Halo 4 


SiO (2-1) 


196.1 ± 0.4 


20.6 ± 0.9 


0.123 


2.69 ±0.11 


3.4 


4.51 ± 0.19 x 10 12 




29 SiO(2-l) 


196* 


21* 




< 0.14 


3.9 


< 0.242 x 10 12 




30 SiO (2-1) 


196* 


21* 




< 0.14 


3.9 


< 0.250 x 10 12 




CS (2-1) 


197.1 ±0.3 


21.9 ±0.6 


0.380 


8.88 ±0.21 


6.5 






CS (3-2) 


196.7 ± 0.2 


22.6 ± 0.4 


0.183 


4.40 ± 0.07 


5.3 






C 34 S (2-1) 


195.2 ± 2.2 


24.8 ± 4.8 


0.022 


0.57 ±0.10 


4.9 


1.96 ± 0.34 x 10 12 




HNCO 


196.9 ±0.2 


19.0 ±0.5 


0.178 


3.59 ± 0.08 


4.7 


106 ± 2 x 10 12 




c 18 o 


206.1 ± 1.2 


7.6 ±4.0 


0.081 


0.66 ±0.18 


10.0 


5.2 ± 2.7 x 10 15 




13 C0 


202.3 ± 0.6 


25.1 ± 1.4 


0.705 


18.8 ±0.9 


25.0 


1.28 ± 0.06 x 10 16 


Halo 5 


SiO (2-1) 


-62.9 ±0.1 


16.3 ±0.3 


0.254 


4.39 ± 0.07 


4.6 


7.37 ± 0.13 x 10 12 




29 SiO(2-l) 


-64.8 ± 0.8 


11.9 ±1.4 


0.032 


0.41 ± 0.05 


3.3 


0.70 ± 0.09 x 10 12 




30 SiO (2-1) 


-61.4 ±2.4 


22.1 ±7.5 


0.917 


0.40 ±0.10 




0.71 ±0.17 x 10 12 




CS (2-1) 


-62.7 ±0.1 


16.7 ±0.1 


0.646 


11.52 ±0.04 


2.9 






CS (3-2) 


-62.7 ± 0.2 


14.4 ± 0.4 


0.528 


8.12 ±0.18 


6.4 






C 34 S (2-1) 


-62.3 ± 0.3 


13.3 ±0.6 


0.074 


1.05 ±0.04 


3.0 


3.62 ± 0.14 x 10 12 




HNCO 


-64.2 ± 0.2 


14.6 ± 0.4 


0.164 


2.54 ± 0.07 


4.3 


75.0 ± 1.9 x 10' 2 


DiskXl-1 


SiO (2-1) 


178.0 ±0.3 


41.3 ±0.7 


0.390 


17.1 ±0.2 


9.1 


28.75 ± 0.39 xlO 12 
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Table A.5. continued. 



Source 


Species 


Velocity Center 


A v 


T A 


J T K dV 


rms 


N 






LSR [km s" 1 ] 


[km s" 1 ] 


[K] 


[K kms- 1 ] 


[mK] 


[cm" 2 ] 




29 SiO(2-l) 


179.4 ± 2.5 


43.4 ±5.0 


0.031 


1.43 ±0.16 


5.4 


2.46 ± 0.27 x 10 12 




30 SiO (2-1) 


188.0 ± 1.0 


24.6 ± 2.6 


0.039 


1.00 ±0.08 


3.3 


1.77 ± 0.13 x 10 12 




CS (2-1) 


176.2 ±0.3 


36.4 ± 0.7 


0.718 


27.81 ± 0.41 


6.3 






CS (3-2) 


174.3 ±0.1 


31.1 ±0.3 


0.591 


19.60 ±0.16 


7.7 






C 34 S (2-1) 


177.5 ± 1.3 


35.1 ±2.8 


0.040 


1.51 ±0.11 


5.1 


5.21 ±0.38 x 10 12 




HNCO 


187.6 ±0.6 


16.4 + 1.8 


0.170 


2.96 ± 0.24 


9.3 


87.4 ± 7.2 x 10 12 




c 18 o 


177.9 ± 6.4 


49.5 ± 13.0 


0.040 


2.12 ±0.46 


13.3 


1.45 ±0.31 x 10 15 




13 CO 


175.6 ± 2.6 


56.4 ±6.1 


0.603 


36.2 ±3.2 


21.2 


2.47 ± 0.22 x 10 16 


DiskX2-l 


SiO (2-1) 


99.7 ± 0.3 


33.2 ±0.8 


0.282 


9.96 ± 0.21 


9.1 


16.72 ±0.36 x 10 12 




29 SiO(2-l) 


95.2 ±3.3 


57.7 ± 19.1 


0.021 


1.26 ±0.25 


5.4 


2.17 ± 0.43 x 10 12 




30 SiO (2-1) 


94.1 ±3.3 


41.4 ±9.8 


0.013 


0.56 ±0.10 


3.3 


0.99 ± 0.17 x 10 12 




CS (2-1) 


98.1 ±0.2 


30.7 ± 0.4 


0.953 


31.10 ±0.37 


6.3 






CS (3-2) 


97.1 ±0.1 


26.0 ± 0.3 


0.580 


16.08 ±0.16 


7.7 






C 34 S (2-1) 


99.7 ± 0.6 


21.0 ±1.8 


0.064 


1.43 ±0.09 


5.1 


4.94 ± 0.33 x 10 12 




HNCO 


96.3 ± 0.5 


30.9 ± 1.0 


0.293 


9.61 ±0.29 


9.3 


283.5 ± 8.6 xlO 12 




c 18 o 


94.5 ± 0.8 


28.9 ± 1.9 


0.194 


5.96 ± 0.34 


13.3 


4.10 ± 0.24 x 10 15 




13 co 


95.3 ± 0.5 


32.1 ± 1.3 


2.056 


70.3 ± 2.5 


21.2 


4.80 ± 0.17 x 10 16 


Disk XI -2 


SiO (2-1) 


69.1 ±0.1 


17.4 ±0.3 


0.522 


9.70 ±0.12 


5.7 


16.28 ±0.20 x 10 12 




29 SiO(2-l) 


68.1 ±2.4 


18.0 ±5.2 


0.020 


0.39 ±0.10 


5.5 


0.67 ± 0.18 x 10 12 




30 SiO (2-1) 


72.5 ± 3.3 


44.5 ± 7.8 


0.018 


0.83 ±0.13 


4.7 


1.47 ± 0.22 x 10 12 




CS (2-1) 


70.0 ± 0.2 


16.1 ±0.5 


1.349 


23.05 ± 0.66 


4.7 






CS (3-2) 


70.6 ± 0.2 


12.9 ± 0.4 


1.320 


18.16 ±0.45 


7.3 






C 34 S (2-1) 


71.1 ±0.4 


14.8 ± 1.0 


0.106 


1.67 ±0.09 


4.7 


5.77 ± 0.32 x 10 12 




HNCO 


68.7 ± 0.4 


18.3 ±0.8 


0.151 


2.94 ±0.12 


7.0 


86.7 ± 3.4 x 10 12 




c 18 o 


68.3 ±2.8 


16.7 ±5.6 


0.058 


1.04 ±0.32 


5.6 


0.71 ±0.22 x 10 15 




13 co 


68.0 ±2.2 


16.6 ±5.4 


0.567 


10.0 ±2.9 


8.0 


0.68 ± 0.20 x 10 16 


DiskX2-2 


SiO (2-1) 


-43.5 ± 2.0 


29.3 ±4.5 


0.037 


1.15 ±0.16 


5.7 


1.92 ± 0.26 x 10 12 




29 SiO(2-l) 


-43* 


29* 


- 


< 0.23 


5.5 


< 0.402 x 10' 2 




30 SiO (2-1) 


-51.8 ±2.8 


28.0 ±4.9 


0.018 


0.54 ±0.10 


4.7 


0.96 ± 0.17 x 10 12 




CS (2-1) 


-43.8 ± 0.3 


19.6 ±0.7 


1.057 


22.09 ±0.71 


4.7 






CS (3-2) 


-43.4 ± 0.4 


21.6 ± 1.0 


0.648 


14.88 ±0.57 


7.3 






C 34 S (2-1) 


-45.0 ± 0.5 


21.0 ± 1.5 


0.093 


2.07 ±0.11 


4.7 


7.15 ± 0.39 x 10 12 




HNCO 


-47.1 ±0.7 


21.7 ± 1.7 


0.087 


2.01 ±0.13 


7.0 


59.2 ± 3.9 x 10 12 




c 18 o 


-46.1 ±0.5 


18.2 ± 1.3 


0.301 


5.81 ±0.35 


5.6 


3.99 ± 0.24 x 10 15 




13 co 


-44.0 ± 0.4 


17.1 ± 1.1 


3.312 


60.2 ± 3.0 


8.0 


4.11 ±0.20 x 10 16 


Diskl 


SiO (2-1) 


54.1 ±0.6 


11.5 ±1.8 


0.040 


0.49 ± 0.07 


3.2 


0.82 ± 0.12 x 10 12 






73.6 ±0.5 


18.0 ± 1.4 


0.062 


1.18 ±0.07 


3.2 


1.98 ± 0.12 x 10 12 




29 SiO(2-l) 


54* 


12 




< 0.06 


2.2 


< 0.105 x 10 12 






74* 


18 




< 0.07 


2.2 


< 0.129 x 10 12 




30 SiO (2-1) 


54* 


12 


_ 


< 0.06 


2.2 


< 0.106 x 10 12 






74* 


18 




< 0.07 


2.2 


< 0.130 x 10 12 




CS (2-1) 


56.8 ±0.3 


14.5 ± 0.6 


0.087 


1.35 ±0.05 


2.2 








74.9 ±0.1 


17.2 ±0.2 


0.327 


5.99 ± 0.07 


2.2 






CS (3-2) 


56.3 ±0.3 


9.8 ±0.8 


0.036 


0.38 ± 0.03 


2.1 








74.8 ±0.1 


14.8 ±0.2 


0.177 


2.78 ± 0.02 


2.1 






C 34 S (2-1) 


50.8 ± 4.5 


12.8 ±7.3 


0.003 


0.05 ± 0.03 


1.9 


0.16 ± 0.10 x 10 12 






73.8 ±0.9 


15.2 ±2.2 


0.016 


0.26 ± 0.03 


1.9 


0.09 ±0.11 x 10 12 




HNCO 


55.3 ±2.8 


16.5 ±6.3 


0.024 


0.43 ±0.15 


4.2 


12.6 ± 4.5 xlO 12 






75.6 ± 0.7 


15.7 ± 1.6 


0.100 


1.66 ±0.16 


4.2 


49.1 ±4.7 x 10 12 




c 18 o 


55.5 ±0.8 


11.5 ± 1.6 


0.075 


0.92 ±0.12 


9.3 


0.63 ± 0.09 x 10 15 






78.6 ± 0.9 


11.6 ±2.3 


0.067 


0.83 ±0.13 


9.3 


0.57 ± 0.09 x 10 15 




13 co 


57.3 ±0.3 


12.6 ±0.8 


0.649 


9.89 ± 0.48 


10.1 


0.59 ± 0.03 x 10 16 






77.6 ± 0.3 


14.4 ± 0.8 


0.647 


9.89 ± 0.48 


10.1 


0.68 ± 0.03 x 10 16 


Disk 2 


SiO (2-1) 


52.3 ±0.7 


18.3 ± 1.8 


0.519 


10.12 ±0.92 


7.9 


17.0 ± 1.5 x 10 12 






76.9 ± 0.3 


19.0 ±0.8 


1.233 


24.94 ±0.91 


7.9 


41.9 ± 1.5 x 10 12 




29 SiOf2-D 


69.0 ± 1.0" 


26.2 ± 2.3 


0.240 


6 67 + 48 


5.0 


11.5 ± 0.8 x 10' 2 




30 SiO (2-1) 


69.2 ± 1.8" 


26.2 ± 4.4 


0.184 


5.13 ±0.73 


5.4 


9.0 ± 1.3 x 10 12 




CS (2-1) 


14.2 ±6.1 


33.0 ±6.1 


0.461 


16.2 ±2.5 


12.9 








55.9 ±6.1 


21.2 ± 6.1 


1.544 


34.9 ± 2.5 


12.9 








79.0 ±6.1 


17.7 ±6.1 


1.922 


36.1 ± 2.5 


12.9 






CS (3-2) 


14.7 ±5.0 


26.9 ± 13.4 


0.317 


9.1 ±3.5 


14.8 








52.7 ± 1.4 


16.4 ±3.5 


1.110 


19.3 ±3.7 


14.8 








76.9 ± 0.9 


17.5 ±2.7 


1.640 


30.6 ± 3.8 


14.8 





C 34 S (2-1) 
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Table A.5. continued. 



Source Species 


Velocity Center 


A v 


T A 


J T K dV 


rms 


N 




LSR[kms-»] 


[km s" 1 ] 


[K] 


[K kms- 1 ] 


[mK] 


[cm" 2 ] 




33.9 ± 2.6 


37.5 ± 5.9 


0.145 


5.79 ± 0.82 


2.9 


20.0 ± 2.8 x 10 12 




71.6 ±0.5 


23.1 ± 1.2 


0.533 


13.12 ±0.75 


2.9 


45.4 ± 2.6 x 10 12 


HNCO 
















67.1 ±0.1 


25.1 ±0.2 


8.868 


236.8 ± 1.5 


11.8 


6982 ± 44 x 10 12 


c 18 o 
















68.0 ±0.8 


23.4 ± 2.0 


0.945 


23.5 ± 1.6 


14.1 


16.2 ± 1.1 x 10 15 


13 co 
















69.1 ±0.6 


28.5 ± 1.5 


5.932 


179.9 ±7.6 


22.2 


12.3 ±0.5 x 10 16 



only possible to fit one velocity components, 
not detected, value for reference taken by the main isotope, 
formal errors of a Gaussian fit. The calibration errors may be larger. 
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Appendix B: Complementary Figures 

This Appendix presents the LVG diagrams of the metastable in- 
version transitions of para-NH3. Most of the sources show two 
kinetic temperature components: one warm which is plotted in 
a red line; and one cool which is plotted in a blue line. In the 
cases where only the warm kinetic temperature component was 
present, the result is showed with a red line. We show the n(H2) 
derived for each component from the CS data, which was used as 
a fixed parameter in the RADEX program. The shaded regions 
correspond to the error associated to each para-NH3 line. When 
a line was not detected we plot their 3 - cr level in a dashed line. 
The error associated to the kinetic temperature and column den- 
sity was estimated computing the x 2 of the line intensities over 
the grid used for the LVG model. We impose Ay 2 = X L ~x\^ m = 1 
which translate in the 68.3% confidence level projected for each 
parameter axis which is showed in a black elipse. 
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Radex, para— NH ;1 Halo 1-1 



Radex, para— NH S Halo 4 




NH,(5,5) 



n(H e )=Q. 50x10* cm" 
Ax 2 = 1 



too 



a oo 



150 
T„, [K] 

Radex, para— NH ; , Halo 1 - 



13 



12 



16 




NH,(5,5) _: 



n(H s ) = 6.71xl0 4 cm" 

AX 2 = 1 



100 



150 
T„„ [K] 

Radex, para— NH a Halo 1-3 




n(E s )= 12.43x10' cm" 
Ax 2 = 1 



150 

T„,» [K] 



Fig. B.l. LVG diagrams of NH3 for each velocity component 
of Halol. Top: 87.7kms _I . Middle: 117.2kms _1 . Bottom: 
138.6 km s" 1 . 




NH.(5,5) 
- NH,(4,4) — ; 



n(H\,)=2.?5xl0 < cm" 
AX 2 = 1 



150 

T» [K] 



Fig. B.2. LVG diagram of NH 3 for Halo 4. 

Radex, para-NH ;! Disk Xl-1 




n(H,,)=12.48xl0 4 cm" 
n(H s )=3. 55x10" cm" 3 
AX 2 = 1 



150 

T» [K] 



Fig.B.3. LVG diagram of NH 3 for Disk Xl-1. The blue lines 
correspond to the low temperature regimen, and the red lines to 
the high temperature regimen. 



Radex, para-NH, Disk X2-1 




n(H ! ) = 5. 62x10' cm" 
n(H 2 ) = 2. 51x10" em- 
Ax 2 = 1 



150 

T» [K] 



Fig.B.4. LVG diagram of NH 3 for DiskX2-l. The blue lines 
correspond to the low temperature regimen, and the red lines to 
the high temperature regimen. 



23 



D. Riquelme et al.: 7\in of X1/X2 orbits interceptions and Giant Magnetic Loops 



18 



Radex, para-NH ;1 Disk XI -2 




n(H 2 )= 12.59x10* cm 3 
n(H,)=5. 62x10* cm" 3 
AX 2 = 1 



150 
T„, [K] 



Fig.B.5. LVG diagram of NH 3 for DiskXl-2. The blue lines 
correspond to the low temperature regimen, and the red lines to 
the high temperature regimen. 
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Fig.B.6. LVG diagram of NH 3 for DiskX2-2. The blue lines 
correspond to the low temperature regimen, and the red lines to 
the high temperature regimen. 
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Fig. B.7. LVG diagrams of NH3 for each velocity component 
of Disk 1. Top: 56.5 km s~'. Bottom: 72.4 km s~'. The blue lines 
correspond to the low temperature regimen, and the red lines to 
the high temperature regimen. 
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Fig. B.8. LVG diagrams of NH3 for each velocity compo- 
nent of Disk2. Top: 48.8kms~'. Middle: 74.5 km s. Bottom: 
95.8 km s . The blue lines correspond to the low temperature 
regimen, and the red lines to the high temperature regimen. 
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Radex CS Halo 1-2 
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Fig. B.10. LVG diagrams of CS for each velocity component of 
«alo 2. Top:-80kms _1 . Bottom: -5 1.9 kms -1 . 



Fig. B.9. LVG diagrams of CS for each velocity component (de- 
fined by the J=2-l CS) of Halo 1 (see Table |A3) . Top: 83.8 
kms -1 . Middle: 117.4 kms -1 . Bottom: 136.9 kms" 1 . The r kin is 
indicated in the lower left corner in each plot. 
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Fig. B. 11. LVG diagrams of CS for each velocity component of 
Halo3. Top: -64.7 km s -1 . Bottom: -13.9kms _1 . 
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Fig. B.12. LVG diagrams of CS for Halo 4. 
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Fig. B. 14. LVG diagrams of CS for each kinetic temperature regime of DiskXl-1. 
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Fig. B.15. LVG diagrams of CS for each kinetic temperature regime of DiskX2- 1 
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Fig. B.16. LVG diagrams of CS for each kinetic temperature regime of Disk XI -2 
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Fig. B.17. LVG diagrams of CS for each kinetic temperature regime of DiskX2-2 



D. Riquelme et al.: T kin of X1/X2 orbits interceptions and Giant Magnetic Loops 



Radex CS Diskl-1 Radex CS Diskl-1 




12 13 14 15 18 12 13 14 15 

Log(N(CS)) [cm"'] Log(N(CS)) [cm"'] 



Radex CS Diskl-2 Radex CS Disk! -2 




12 13 14 15 16 12 13 14 15 

Log(N(CS)) [cm-"] Log(N(CS)) [cm-'] 



Fig. B. 18. LVG diagrams of CS for each kinetic temperature regime, for each velocity component of Diskl. Top: 56.8 km s 
Bottom: 74.9 km s~'. 
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Fig. B. 19. LVG diagrams of CS for each kinetic temperature regime, for each velocity component of Disk2. Top: 14.2 km s 
Middle: 55.9 km s -1 . Bottom: 79.0 km s -1 . 



